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PREFACE 


The power and adaptability of the matrix as a mathematical tool is 
evident by matrix techniques becoming used in an increasing number of 
fields. No mathematical education is now complete without some 
knowledge of their use. 

The material in this book has come from a course in matrices and 
their applications given to students training as mathematics teachers. In 
each section, the starting point is an application or situation in which 
matrices arise. Techniques are then developed for dealing with the 
situation, and exercises and investigations are provided through which 
the reader may extend his understanding of the ideas. 

We hope that students in sixth forms, Colleges of Education and of 
Technology, and in the first year at University will find this material 
useful. Teachers will find material suitable for extending their own 
knowledge of the subject, and also some situations which could form the 
basis of lessons for secondary school children. 

We are indebted to our own students on whom we have tried out this 
material; to Dr. T. J. Fletcher (H.M.I.) for many of the ideas on which 
Chapters 7 and 8 are based; to the editor of the Monograph series, 
Mr. A. J. Moakes, for his generous help and encouragement and to the 
publishers. 

Nottingham College of J. R. Branfield 

Education, 1969. A. W. Bell 
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MATRICES AND THEIR ALGEBRA 

1.1 Introduction 

Matrices were invented in the context of sets of linear equations, a 
topic which we shall discuss later, but the range of their applications has 
increased enormously in recent years, particularly in the planning of 
industrial operations, of production and marketing. We shall introduce 
matrices in relation to some of these newer applications. 


1.2 The Matrix 

There are many situations where it is natural to present data in the 
form of a table. One of these is shown in Table 1.1. 



English 

Maths 

French 

Men 

45 

50 

15 

Women 

55 

20 

20 


Table 1.1. The matrix F. First-year students taking English, 
Mathematics or French at Nottingham College of Education. 

A matrix is any rectangular array of numbers. The array of numbers 
obtained by omitting the headings in Table 1.1 is a matrix of two rows 
and three columns; this is called a ‘2 x 3’ (two-by-three) matrix and will 
be denoted by F; we shall write it with round brackets thus: 

/45 50 15\ 

\55 20 20y* 

A matrix consisting of only one row or column is called a vector ; thus 

/45\ 

the column ( 55 ) would be called a column vector, and the row 

(45, 50, 15) would be called a row-vector. The individual numbers of a 
matrix are called its entries (or components). 

1 
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1.3 Vector Addition 

It is the object of matrix theory to develop methods of handling 
matrices as single units instead of having to deal with all their entries 
separately, in particular to give suitable meanings to addition and 
multiplication for matrices. This may be compared with what we do 
with fractions; the fraction £ consists of the pair of integers 3, 4; but we 
have developed ways of dealing with a fraction as a single unit, so that we 
know what we mean by the sum of two fractions, we know when to use 
this sum, and we know also how to obtain the sum from the two pairs of 
integers which constitute the original fractions: (3, 4)+ (1, 2) = (5, 4). 
Just as a fraction is thus a more complicated but more widely useful 
type of number than an integer, so a matrix is a more complicated but 
more widely useful type of number than either. We shall now consider 
what mathematical operations might be performed on the matrix F, and 
hope to see whether any of these might usefully be regarded as addition 
or multiplication for matrices. 

What further information can we gain from F ? Can we find the total 
number of students involved? Yes, this can be done by adding the 
numbers of men and women in each column, to give column totals; 
English 100, Mathematics 70, French 35, and grand total 205. This 
result may be checked by finding the row-totals first: 

men 45+50 + 15 = 110, 

women 55 + 20 + 20 = 95, 
grand total 110 + 95= 205. 

These totals are included in the enlarged matrix F 2 shown in Table 1.2. 

English Maths French Row Totals 




e 

a 

f 

r 

Men 

m 

45 

50 

15 

110 

Women 

w 

55 

20 

20 

95 

Column Totals 

c 

100 

70 

35 

205 


Table 1.2. Table 1.1, with totals. The matrix F 2 . 


Let us denote the row and column vectors of this matrix as follows: 

Rows: Men m, Women w, row totals r 
Columns: English e, Maths a, French f, column totals c 

We shall define the sum of two vectors as the vector obtained by 
adding corresponding entries of the two vectors, and we can then write 
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what we have done as m+w = c, this being a brief statement for 
(45, 50, 15, 110)+ (55, 20, 20, 95) =(100, 70, 35, 205) and e + a+f=r 

which means ( 5s) + (V) + (V) = ( 95V 
VlOO/ \70/ \35/ \205/ 


Exercise 

1.1 Compute the following vectors: 

e + a, e-a, r-e, m + c, m - w. 

The information given in Table 1.2 is lacking in one important respect 
- quite apart from the fact that it only quotes figures for three subjects. 
At this college, some students take two subjects, while others take only 
one. The row-totals do not therefore represent the total numbers of 
students involved in the three subjects, since those who take, for 
example, English and Mathematics will have been included twice, once 
under each heading. Table 1.3 shows the numbers of these, and omits 
the French students, for simplicity. 

both English 


English Maths 

and Maths 

Row ‘totals’ 

e a 

b 

*3 


Men 

m 3 

45 

50 

2 

93 

Women 

w 3 

55 

20 

3 

72 

Column Totals 

C 3 

100 

70 

5 

165 


Table 1.3. The matrix F 3 . 


Here the row totals are calculated so as to give the total number of 
students involved; the number taking both subjects must be subtracted 
from the sum of the separate English and Maths lists, since they will 
have been counted twice. In vector notation, we have 

r 3 = e + a - b 

using the letters indicated in the table. 


1.4 Matrix addition 

Table 1.4 shows the numbers of second-year students in the same 
categories as those of Table 1.1 (this was for first year students). The 
matrix obtained by omitting the headings here will be called S, 
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English 

Maths 

French 

e 

a 

f 


Men m 

25 

30 

0 

Women w 

35 

18 

0 


Table 1.4. The Matrix S. 


We shall define the sum F + S of the two matrices F and S as the 
matrix formed by adding entries in corresponding positions, thus: 

/45 50 15\ (25 30 0\ /70 80 15\ 

\55 20 20/ + \35 18 0/ \90 38 20/* 

This gives us a matrix showing the total numbers in the first and 
second years combined, in each of the six categories. 

1.5 Multiplication of a Matrix by a Scalar 

Following a natural development we write F + F as 2F, F + F + F to n 
terms as nF ; and we extend this definition to apply to multiplication by 
any number k drawn from the same set as the entries of the matrix thus: 
The matrix kA is that obtained from A by multiplying each entry by 

k. Thus 

0 0*2\ _/ 0 0-62\ 

-1 3 ) 9-3 )■ 

This is called multiplying a matrix by a scalar; a single number is called 
a scalar to distinguish it from a vector and a matrix. 

Exercises 

Compute the following matrices: 

l. 2 S 3 (S modified in the same way as F was to give F 3 ); the b 

column is 

1.3 F - S, F 3 - S 3 , 4S. 

1.6 Matrix multiplication 

This refers to the multiplication of one matrix by another. We shall 
begin by investigating whether there is any way of combining two 
vectors which might reasonably be called multiplication. 

Suppose the first-year students of Table 1.2 are to be entered for 
examinations. The examination fees per student are 40p for English, 
30p for Mathematics (because it’s easy to mark!) and 50p for French (to 
cover the extra cost of an oral). What will the total cost be? Using the 
bottom row of the table, it is, in pence, 
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100.40 + 70.30 + 35.50 = 4000 + 2100 +1750 = 7850 (£78*50). 

This is one possible way of multiplying the vectors (100, 70, 35), (the 
subject-numbers vectors) and (40, 30, 50), (the fee vector); it can be 
described as multiplying corresponding entries and adding the results. 
This product of two vectors is called the inner product, or sometimes the 
scalar product , since the result is a single number or scalar. The inner 
product of vectors a and b is denoted by a . b, which accounts for yet a 
third name, namely the dot product of the vectors a and b. 

Exercises 

1.4 If all students have to pay a games subscription, £3 per year for 
men and £2 per year for women, what total sum is paid by (a) 
English students, (b) Maths students, (e) all these students? (Use 
Table 1.3, F 3 .) Show how these results may be obtained as inner 
products of the subscription vector s = (3, 2) with the appropriate 
vectors taken from Table F 3 . 

1.5 If a = (2,1, 0), b = (- 3, 0, 2), c = (0,1, 0), form the inner products 
a . b, b . c, a . c. 

Continuing the discussion of the last paragraph, we might wish to 
calculate the examination fees to be paid by men and women separately 
as well as the total. We should then be forming the inner product of the 
fee vector v = (40, 30, 50) with each of the rows of the matrix F 2 (Table 
1.2), omitting the row totals, and recording the three results. 

We have, for the men’s fees, 

(45, 50,15). (40, 30, 50) =45.40+ 50.30 +15.50 
= 1800 + 1500 + 750 
=4050 

and for the women’s fees 

(55, 20, 20). (40, 30, 50) = 55.40+ 20.30+ 20.50 
=2200 + 600 + 1000 
= 3800 

and a total of 7850, as before. We record the results as a column vector 
/4050\ 

( 3800 ] , and the whole calculation would be written in matrix notation 
'7850' 

/ 45 50 15\ /40\ /4050\ 

as ( 55 20 20 (30 = 3800 . 

M00 70 35/ '50/ '7850' 

/4050\ 

Here we may see good reason for writing the vector I 3800 I as a 

M850/ 



6 


MATRICES AND THEIR APPLICATIONS 


column, since the components refer to men, women and total, and these 
are the titles of the rows of the matrix F 2 ; but it is not obvious why the 
fee-vector (40, 30, 50), whose entries refer to the three subjects English, 
Maths, French should be written as a column rather than a row. The 
real answer to this cannot be given until we have more experience of 
matrix multiplication, but at least it can be seen that there must be 
some agreement about whether the vector (40, 30, 50), when written 
following a matrix, is to be multiplied by the rows or the columns of the 
matrix. The agreement is that multiplication is always row by column, 
that is, each row of the first matrix is multiplied by each column of the 
second. (In this case, the ‘second matrix’ is a single column vector.) Also, 
in the product matrix, the row position of each product entry is that of 
the first matrix from which it was derived and its column position is 
that of the column of the second matrix from which it came. This may 

f w \ 

be seen more clearly if we put another column alongside I 30 1 to make a 

'50' 

3x2 matrix, and perform the multiplication: 

/ 45 50 15\ /40 1\ /4050 -5\ 

55 20 20 ( 30 -1 =1 3800 35 . 

M00 70 35' '50 O' '7850 30' 


Here the second column of the product matrix is obtained from the 
second column of the second matrix, and the row position of the entries 
in that product column is determined by the row of the first matrix 
from which they came. 

In practice, we work like this: 


then 


then 


/45 50 15\ / 40 *\ 45.40 

* * * | I 30 *U+ 50.30 

\ 50 */ +15.50 



745 50 15\ 


* * * 

* * * / 


* 

* 

* 



45.1 

+ 50.-1 
+ 15.0 


'4050 




40 

30 

50 


*\ 55.40 

*)-> + 20.30 
J +20.50 




and so on. 

We define A 2 as A . A and say that two matrices are equal if their 
corresponding entries are equal. 
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Exercises 

1.6 IfA = (j _f),B = (J 2) ’ C = (o l)’ D = (l o)> compute 

the products of all possible pairs of these matrices, and A 2 , B 2 , C 2 . 

1.7 If 



/I 0 0\ (0 1 0\ 

G = ( 0 1 O), H=(l 0 0 , 

Vo 0 1/ Vq 0 V 

compute the products of all possible pairs of these matrices, and E 2 , F 2 . 

Note the effect of C, D, G, H on the matrices they multiply. 

1.7 Other ways of combining matrices 

In this chapter so far we have taken a table of numbers and have 
discussed what mathematical operations might be performed in it, with 
a view to showing the operations of matrix addition and multiplication. 
We now look at another table with a view to seeing what mathematical 
operations are suggested - whether matrix addition and multiplication 
yield anything of interest, and whether any other operations emerge. 

Table 1.5, taken from Volume II of the Crowther Report, shows the 
number of candidates obtaining scores 4, 3, 2, 1, 0 in an ability test, 
classified according to the number of children in the family in which 
they were brought up. The purpose of this table, in its context, was to 
demonstrate that the mean score showed a steady decrease as the size of 
family increased; the last row is the relevant one. In the calculation of 
this row a number of matrix operations can be observed. The first 
figure in the total score row is 4.5 +3.14 + 2.9 + 1.24 + 0.13; i.e. c . c ly 
and if the whole inner 5x6 matrix, omitting the headings and totals, is 
denoted by M, the row t 3 is the matrix product cM. (Verify this.) The 
row t c can also be written as a matrix product, if desired, as uM, where 
u is the row vector (1, 1, 1, 1, 1), and this might be useful in enabling 
both calculations to be dealt with by the same machine programme. 
The derivation of m from t c and t s , by component-wise division, does 
not correspond to any matrix operation so far defined; we are obtaining 
a vector (c v c 2f ... ,c n ) from (a ly a 2i ... y a n ) and (b ly b 2> ... y b n ) by putting 
Ck = a klbjcf°* all A 

Expressing this in terms of multiplication, we have a product bvc say 
(we invent the symbol v, calling it vee), defined for any two vectors of 
the same number of components by (bvc) k = b k c k . This easily generalises 
to matrices of the same number of rows and columns and gives us a 
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c 

Score in 
ability test 

Cl 

1 

C 2 

I' 

2 

C 3 

s[o. Of d 

3 

I b C 4 
lildren: 

4 

C 5 

in famil 

5 

c 6 

y 

6 or 

more 

t 

Totals 

s 4 4 

5 

4 

4 

2 

i 

i 

17 

s 3 3 

14 

24 

32 

23 

16 

22 

131 

s 2 2 

9 

23 

31 

23 

15 

27 

128 

Si 1 

24 

31 

34 

39 

37 

42 

207 

s 0 0 

13 

26 

39 

40 

36 

76 

230 

t c Total cands. 

65 

108 

140 

127 

105 

168 

713 

t s Total score 

104 

165 

208 

162 

119 

166 

924 

m Mean score 

1-60 

1*53 

1-49 

1-28 

M3 

0-99 

1-30 


Table 1.5. Number of candidates obtaining different scores in 
an ability test, classified by family size. 

(From the Crowther Report, Vol. II, p. 178.) 


matrix product, MvN say, which corresponds to the matrix sum M + N 
which we have already defined, multiplication simply replacing addition. 
The question of whether this new element-wise product MvN is more 
or less valid than the previous row-by-column product MN we leave 
over until we have completed our examination of the table above. 

Exercises 

1.8 Consider how each figure in the final totals column t is obtained, 
noting what matrix operations are used. 


1.9 Make a similar study of Table 1.6, calculating totals and observing 
what operations are involved. 


No. of 
children 

Professional 

and 

managerial 

Clerical 

Skilled 

Semi¬ 

skilled 

Unskilled 

Total 

1 

136 

120 

403 

86 

65 


2 

277 

201 

772 

172 

108 


3 

150 

155 

621 

176 

140 


4 

78 

77 

442 

108 

127 


5 

20 

43 

261 

98 

105 


6 or more 

31 

60 

424 

163 

168 



Table 1.6. Number of families of different sizes in different 
socio-economic groups. 

(From the Crowther Report, Vol. II.) 
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1.8 Laws of matrix Algebra 

We now have one way of adding and two ways of multiplying 
matrices, so that we have made some progress towards our aim of 
enabling matrices to be manipulated like numbers. It is also clear from 
the exercises on p. 7 that the laws of matrix algebra differ from those of 
numbers in some respects - two matrices can only be combined if they 
have the right numbers of rows and columns, and even if matrices A 
and B can be multiplied (by matrix multiplication) in both orders, AB 
and BA, the product matrices are not in general the same. 

At this point we must decide which of the two alternative multiplica¬ 
tions to adopt, and we choose the row-by-column one as this is the one 
which is relevant to the applications we shall consider. 

We can now summarise the laws of this matrix algebra. 

First, it is convenient to define the multiplication of a matrix by a 
number: if A is a number belonging to the same number system as that 
from which the entries of the matrix M are drawn, we define AM to be 
the matrix obtained from M by multiplying every entry by A. 

The laws are as follows: 

Associative Laws for Addition and Multiplication 

A + (B + C) = (A + B) + C; A(BC)=(AB)C. 

Commutative Law (addition only) 

A+B=B+A 

Distributive Law, and laws concerning multiplication by a number: 

A(B + C) =AB + AC 

A(A + B)=AA + AB, (A+/x)A=AA +/xA 

A(AB) =AA . B, (Aju)A=A(/xA). 


All of these can be verified from the definitions, most of them easily. 
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TRANSFORMATIONS OF THE PLANE 


2.1 Use of 2 x 2 matrices 

Matrices have an important part to play in the study of geometrical 
transformations. In this chapter we consider various transformations of 
the plane using 2x2 matrices; the investigation leads us to a number of 
important matrix ideas. 

Consider a square piece of elastic material OABC held at the origin of 
coordinates O (Fig. 2.1). It has a sleeve along OA and OC (through 
which the axes pass) so that when the material is pulled in the x and y 
axis directions, the material can move along the axes but will stay 
attached to them. The material is stretched by a factor of 2 in the x 
axis direction and by a factor of 3 in they axis direction. 



By stretching, the square OABC is transformed into the rectangle 
OA*B*C*. 

A( 1, 0)^>A*(2, 0) 

B( 1, l)-^*(2, 3) 

C(0, 1)->C*(0, 3) 

and O remains fixed. 


In general, a point P(x, y) in the material will be transformed into the 
point P*(x* y y*) where 


x* = 2af) 

y* = 3yj 

10 


• (i) 
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Equation (1) can be written as 

x* = 2x + 0y 
y* =0x + 3y. 

These may be put into matrix form as 



- which may be verified by performing the matrix multiplication on the 
right-hand side. 

Thus we see a close relationship between the matrix ^ ^ and the 

transformation of the plane by stretches of factors 2 and 3 in the x and y 
directions respectively. In this chapter we shall be exploring this 
relationship, considering different transformations and the matrices to 
which they correspond. We shall be able to use matrix algebra in 
geometry whenever we have a transformation which can be considered 
as a set of linear equations. Such transformations of the plane as 
reflection in a line, rotation about a point, translation and so on can be 
treated as a set of linear equations and can thus be put into matrix form. 
Matrix algebra can be interpreted in geometrical terms whenever such 
interpretation is useful and helpful, and similarly geometrical trans¬ 
formations can be treated and manipulated by matrix algebra. 

2.2 A matrix to rotate points through 90° 

A point P(pc , y) is rotated through 90° anticlockwise about the origin 
O to the point P*(x *, y*). The point P is mapped onto the point P # by a 
transformation of the plane which rotates points through 90° about the 
origin. From the geometry of the Fig. 2.3 


x * = -y 



Fig. 2.3 
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The matrix form of these equations is 



Hence the matrix which will rotate points of the plane through 90° anti¬ 
clockwise about the origin is 

G -*)• 



and the 


For example the point becomes ^ ^ ^ 

point ^ becomes ^ "" o) ( ~ l) = ( -1)‘ matrix (i ~~ o) 

has the property of rotating all points of the plane through 90° about O, 
keeping the coordinate axes fixed, and can be considered as a transfor¬ 
mation which rotates the plane through 90°. 


2.3 A matrix to rotate the plane through any angle a 

The plane is to be rotated through an angle a (anticlockwise) about 
the origin O, so that P(x,y ) goes to P # (# # , y # ) with OP = OP* and 
lPOP* = a. What is the matrix to perform this transformation? 

To find a matrix which will perform a particular transformation we 
need to express the coordinates of the new position of a general point of 
the plane, in terms of its coordinates before the transformation; this 
requires a pair of equations. 

Let the angle which OP makes with the x axis be 9 ; and let OP and 
OP * each be of length r (Fig. 2.5). 

Then x*=r cos (a + 9) 

=r cos a cos 9 -r sin a sin 9 ; 
which since x = r cos 9 and y = r sin 9 
becomes x* =x cos a -y sin a. 



TRANSFORMATIONS OF THE PLANE 


13 


Similarly y * = r sin (a + 0) 

= r sin a cos Q + r cos a sin 0 
and so y* =x sin a+y cos a. 

Hence the pair of equations is 

x* = x cos a -y sin a 
y* =x sin a +y cos a 
from which we derive the matrix form 

( x*\ _ /cos a - sin a\ / 

\sin a cos a/ \y/ * 

Thus the geometric transformation represented by the matrix 

cos a - sin a 
sin a cos a 

maps the point P(x , y) onto the point y # ) by a rotation of the 

plane through an angle a about the origin. Putting a = 30° gives the 
matrix to rotate the plane anticlockwise through an angle of 30° about O 





Considering the transformation by this matrix of a rectangle with 
vertices 0(0, 0); A( 2, 0); B( 2, 1); C(0, 1), we find the positions of the 
vertices of the transformed figure by using 
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Putting in (2, 0) for (, x , y) we obtain the coordinates of A * as 



that is x*=J3 and y* = 1, so A( 2, 0)->A*(J3 9 1). The reader should 
plot these coordinates and verify that the rectangle has been rotated 
anticlockwise through an angle 30°. 

An alternative way of deriving the matrix for a rotation about the 
origin is as follows: 


Consider the matrix 


“(* $)• 

This maps into (^j , and into . 


But under the rotation, (see Fig. 2.6) 


(lW COSa ) and (?W" Sin “V 

\0/ \sma/ \1/ \ cos a/ 

The only possible matrix which can produce this transformation is 
therefore 

( cos a — sin a\ 
sin a cos a/ * 

Moreover, this matrix rotates every point of the plane about O in the 
same way, since 


( cos a 
sin a 


-sin a 
cos a 


(x\ _/cos oc 
\jy/ \sin a 


= x 


( cos a 
sin a 


-sin 

cos 


-MM?) 

ina\/l\ /cos a 
os a/ \0/ + ^\sin a 


-sin 

cos 


- sin a\ /0\ 
cos a/ \1/ 


using the laws of matrix algebra (p. 9) 


= xf C P S a> ) +y( S ^ n 
\sin a/ y \ cos a/ 

so that the right-angled triangle ONP has rotated into the position 
ON*P *, and thus OP has turned through the angle a about O. (Fig. 2.7). 
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2.4 Movement of figure or movement of axes ? 


The equation 

A/3 

i\ 



( 3 - ? 

2' 

JZj 

© 


' 2 

2' 



describes a transformation of (x, y) into (x*,y*); for example, if 
(«xr, y) = (2, 0), then (x*> y*)=(J3 y 1). In the section above we have 
regarded (#, y) as denoting points of the rectangle in its initial position, 
and (x*> y*)> as the coordinates of the same points after the rectangle 
has been rotated. We might, however, adopt a different point of view 
and regard (x* y y*)> as the coordinates of the points of an unmoved 
rectangle with respect to new axes\ axes which are obtained by rotating 
the old ones through 30° clockwise . The same matrix equation may 
represent either a rotation of the rectangle through +30°, the axes 
remaining fixed, or a rotation of the axes through - 30°, the rectangle 
remaining fixed. To put this another way, the coordinates ( x,y) 
describe the position of points of a figure with respect to a given pair of 
axes. If (#, y) changes, then the position of'the figure with respect to 
the axes changes; this may be regarded as a movement of the figure, or 
of the axes. In what follows a movement ‘of the plane’ means a movement 
of the figures in the plane, the axes remaining fixed. Thus the matrix 
for a rotation of the axes through a positive angle a is obtained from the 

equations on page 13 by changing a to - a; this gives ^ ^ ^ * 

as the matrix which keeps points of the plane fixed and refers them to 
new axes rotated through an angle a anticlockwise about O. 


2.5 Other transformations in the plane 

At the beginning of this chapter we considered the matrix 

which may be described geometrically as a double-stretching of the 
plane, since for each point its distance from the y axis is multiplied by 
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2 and its distance from the x axis is multiplied by 3. This type of trans¬ 
formation will be called a ‘double stretch’. The matrix . . 

\sin a cos a 

has the geometric property of rotating the plane through an angle a 

about O and the matrix ( C ? S a Sln a ) describes the plane with 
\-sin« cos a/ r 

respect to axes which are rotated through an angle a about O from an 
original set of axes. Other transformations of the plane to be considered 
are reflection, enlargement, shear and projection. The first two of these 
are straightforward and are dealt with in the exercises following. 


Exercises 

2.1 By finding suitable linear equations, find the matrix which will 
transform the plane so that each point of the plane 

(a) has its distance multiplied by a factor a ; y ordinate remains 
constant (linear stretch); 

(b) has its distance from the origin multiplied by a factor k 
(enlargement); 

(c) is reflected in they axis (reflection). 


2.2 Investigate the transformation of the rectangle OABC with 
vertices A( 2, 0), B(2, 1), C(0, 1) when the plane is transformed by 
the matrices 



2.6 The shear transformation 

A shear may be considered as a transformation of the plane in which 
each point is moved along a line parallel to the x axis a distance propor¬ 
tional to its y-value. For the point P(x,y)> if a is the proportionality 
factor, then P is moved parallel to the x axis a distance ay . Thus 
x*=x + ay> and since y* =y, the matrix which describes this shear 

transformation is 

Graph paper is needed in the following exercises in this chapter - it is 
important that the transformation of the plane by matrices is performed 
and seen. 
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Exercises 

2.3 Using the shear matrix ^ ^ show the transformation of the 

figures OABC, ABDE and BDFG, where O is the origin, A(2, 0); 
B( 2, 1);C(0, l);D(-2, l);E(-2, 0);F(-2, -l)andG(2, -1). 

2.4 Find the matrix which transforms the plane such that each point 
is moved along a line parallel to they axis by a distance proportional 
to its x- value. 

2.5 Investigate the transformations of the plane by the matrices 

(i) and (ii) ^ ^ - consider the transformation of the 

figures in 2.3. 

2.7 Projection 

A ‘projection' will be considered as a transformation of the plane in 
which points are projected parallel to an axis onto a line through the 
origin. 



For the point P(x y y) its x coordinate remains unchanged in a projec¬ 
tion parallel to the y axis. The y coordinate of P m becomes mx for 
projection onto the liney = mx and so 

x* =x 
y # = mx, 

and the matrix which describes a projection of the plane onto the line 

y = mxi& ($n o)- 

2.8 Unit and zero matrices 

These two important matrices are connected with (i) the transforma¬ 
tion which maps each point onto itself and (ii) the transformation which 
maps all points into the origin. 
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(i) If every point maps onto itself we have 


and 


x* =x = l . # + 0 .y 

y*=y = 0 . # + l . y , 


so that the matrix corresponding to this transformation is 



This is called the unit matrix and is denoted by I. In matrix algebra 
it plays the same part as the number 1 does in ordinary algebra. The 
mapping it performs is called the identity mapping. 


(ii) If every point P(x , y) is mapped into the origin (0, 0), we have 


and 


# # = 0 . # + 0 .y 
y* = 0 . x + 0 . y. 


Thus the matrix which has the property of mapping all points into the 
origin is ; it is called the zero matrix. 


2.9 Finding the transformation of a given matrix 

Instead of finding a matrix to perform a particular transformation, we 
can also find what transformation a particular matrix performs on points 
of the plane. 


Consider the matrix 



in particular, let us see how it trans¬ 


forms the rectangle OABC. (We could use a square, but this may be too 
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regular to show the transformation - so, of course, may be the rect¬ 
angle.) The equation for the transformation is ^ 0 l) (j;) ’ so 

that the point B( 2, 1) will be given by = ^ 0 l) (l) * t ^ iat * s 

x* = - 2 and y* = 1, so that B( 2, 1) is mapped onto #*( - 2, 1). Similarly 
A(2, 0)->^4*( - 2, 0) and C(0, 1)->C*(1, 0), Fig. 2.10. The general 

point P(x, y) becomes P*(x , -y), so the matrix ( \ has the 


geometric property of reflecting points of the plane in the axis of y. 


2.10 Summary of plane transformations and corresponding 2x2 
matrices 


1. identity transformation 


(Unit matrix I) 


2. zero transformation ^ (Zero matrix) 

3. reflection (in x axis) ^ 

4. rotation of figures through a about O 

r> . . • , r \ ( cos a sin a\ 

5. rotation (or axes) ( 

v 'V^-sina cos a/ 

6. enlargement (O fixed) ^ 

7. linear stretch (parallel to x axis) ^ ^ 

8. double stretch ^ 

9. shear (parallel to x axis) ^ ^ 

10. projection (parallel toy axis, onto liney = mx), 


a — sin a 


m 0/ ’ 


Exercises: 

(Graph paper needed.) 

2.6 (a) State the matrix which will enlarge a figure k times. Show 

that this matrix is k times the unit matrix. 
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(b) Show diagramatically what happens to the unit square with 
vertices (0, 0); (1, 0); (1, 1); (0, 1), when k = 2,\, -3. 

(c) Instead of the unit square find what transformation takes 
place, when k =2, when the figures are 

(i) rectangle (0, 0); (2, 0); (2, 1); (0, 1) 

(ii) circle x 2 +y 2 = 1 

(iii) square(-1, — 1); (1, -1); (1, 1); (-1, 1). 

2.7 For the matrices 

c . ?) (-0 ?) (i -?) (-1 -?) 

n (? (-{o) (-? -j) 

show their effect on the rectangle in 2.6(c), (i) and identify the 
transformation of each of the matrices as one operation. 

2.8 For the linear stretch matrix ^ ^ let a = 2, - 3 and show the 

effect on the unit square and figures in 2.6(c), (i) and (ii). 

In the instance of the circle, prove algebraically that the trans¬ 
formed circle is an ellipse. 

2.9 Investigate the transformations described by the following 
matrices 



2.11 Transformations combined 

In this section we consider whether ways of combining matrices - 
by addition, matrix multiplication and element-wise multiplication - 
give rise to any corresponding combination of the transformations which 
the matrices describe. It turns out that only one mode of combination is 
significant here - matrix multiplication. 

Consider a point P(x , y) rotated through 45° about O to P*(x* t j>*) 
and then reflected in the axis of x to P**(#**, y**), Fig. 2.11. From P to 
P* we have (from the rotation of the plane matrix with a = 45°) 
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From P* to P** we have 

1 0\/x*\ 

\y**J \o -lAW' 


• ( 2 ) 


We seek now the single matrix which will describe the whole transfor¬ 
mation. 


Substituting from (1) for 



into (2) then 




( 3 ) 


/1 

_1 \ 

_ V2 

A 

A- 1 

- 1 / 

V2 


1 ( 1 

-1 

II 

1 

H-* 

-1 


0 (using the row on column method 
for multiplying the matrices) 



Thus the matrix which will perform a rotation through 45° followed 
by a reflection in the axis of x is 


J2 



-1 

-1 


)• 


Checking on the point P( 1, 1), Fig. 2.12, a rotation of 45° followed by a 
reflection in the axis gives P** as (0, — J2). From the matrix we get 


V 1 

V2\-l 
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The general proof that row on column matrix multiplication is needed 
is considered in the next section. 


Exercise 

2.10 Find the matrix corresponding to the performance of the two 
transformations above, in the reverse order. 


2.12 Matrix multiplication and the combination of transforma¬ 
tions 


We see from the above exercise that the order of the matrices in (3) 
is important. Moreover, if A is the rotation matrix and B the reflection 
matrix, then from (3) the combined effect of performing A followed by 
B is given by BA; the matrix applied first being written second. (This 


has come about by the substitution used in placing 



from (1) into 


(2).) For three transformations described by matrices A, B and C the 
combined effect of the matrices A followed by B followed by C will be 
given by computing CBA. 

The process of matrix multiplication is directly related to following 
one transformation by another and is sometimes defined from this 
point of view. The working is as follows. 


Let A be the matrix 



which transforms the point P(xy) into 


P*(x*> y*). The transformation can be written as a set of linear equa¬ 
tions, 


x* = ax + by 
y* = cx + dy 


• ( 4 ) 


in matrix form 



or as x* = Ax. 
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The transformation represented by A is followed by a transformation 


B, where B is the matrix . By the matrix B, the point P*(x*, y*) 

is mapped onto P**(x** t y **), and this can be written in matrix form as 
x**=Bx* 


From x* = Ax and x # * =Bx* then if we wish to obtain x** in terms 
of the original coordinates x, then substitution takes place and 
x** = Bx* =B(Ax) = (BA)x. Note that the associative law for matrix 
multiplication is used here. 


Now the transformation B can be put into a set of linear equations as 


x** =px* +qy*\ ,j-\ 

y** = rx* + sy*j . 

Using equation (4) to find x** andy** in terms of x andy 


x** =p(ax + by) + q(cx + dy) 
y* * = r(ax + by) + s(cx + dy), 

that is 

x* * = (pa + qc)x + (pb + qd)y\ 
y ** = (ra + sc)x + (rb + sd)y J 


■ ( 6 ) 


The combined effect of the matrix A followed by the matrix B is 
given by equation (6). These in matrix form are: 


/***\ 

Jpa + qc 

pb + qd\ / x\ 

vv**/ = 

~\ra+sc 

rb+sd) \y) 

From this result and x** = 

(BA)x 


/ x 

\j;** 

II 

■G* i)0 


we have 


(P 4\( a b\_/pa + qc pb + qd\ 

\r sj\c d) \ra + sc rb + sdj 

giving the row on column method of matrix multiplication. 


Exercises 

(Graph paper needed.) 

2.11 If M.(-J ?) Q-(J “J) 

evaluate the products QM and MQ. What transformations do these 
four matrices represent? 
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2.12 (a) With the matrix which will rotate points of the plane through 
30°, show the transformation of the unit square. Transform this 
figure by a matrix which shears by a factor of 3 in the x axis 
direction. Find the matrix which will rotate points through 30° 
followed by a shear of factor 3 in the x axis direction. Use this 
matrix to transform the unit square. Show diagramatically that the 
result of the separate transformations and of the combined trans¬ 
formation is the same. 

(b) Find the matrix which will rotate points through 45° followed 
by a reflection in the x axis. Find the matrix which will reflect 
points in the x axis and then rotate them through 45°. Are the 
resultant matrices for these two sets of operations the same? 
Explain diagramatically. 

(c) Find the matrix which will enlarge the lengths of a figure to 
twice their size, then stretch it in the direction of the y axis by a 
factor of 3 and then reflect it in the line x=y. 

2.13 Find the single transformation which is equivalent to reflection 
in x =y followed by reflection in the x axis; 

(a) by sketching what happens to a figure. 

(b) by multiplying the corresponding matrices. 

2.14 Find the matrix which will rotate the plane first through an angle 
a and then through an angle j8. The resultant matrix will be 
equivalent to the matrix which rotates the plane through an angle 
a + j8. Hence obtain the ‘addition formulae’: 

cos (a 4- £) = cos a cos £ - sin a sin /? 
sin (a 4- j8) = sin a cos j8 4- cos a sin 

2.15 Find the matrix A which rotates points of the plane through 90°, 
and the matrix B which rotates points of the plane through 180°. 
Compute AB and BA and comment on your results algebraically 
and geometrically. Find the matrix M which reflects points of the 
plane in the line y = x and the matrix N which reflects points of the 
plane in the y axis. Compute MN and NM and comment on your 
results algebraically and geometrically. Is there any relationship 
between the matrices AB, MN, NM - if so, why? 

2.16 If A = ^ 2)'’ B = (l 0)’ C = (o -l) ; investi S ate and 

comment on the transformations represented by AB, BA, A(BC) 
and (AB)C. 
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2.17 (a) Show the effect of the matrix ^ ^ m transforming the 

unit square. Can the matrix be put in terms of the combined effect 
of basic matrices ? 

(b) Investigate the transformation effect of the matrices 

(? g i> (-j:) 

(Note how the vertices correspond). 

2.18 The reflections A, B take the point (, x , y) into (;y , x) and (#, -y) 
respectively. Write down the matrix for each transformation. 

Describe geometrically the effects of AB, BA, AB 2 , (BA) 2 , ABA 
and BAB. 

Verify that 

(AB) 4 = I = (BA) 4 = (ABA) 2 . 

Write ABABA and BABABA using in each case four or less of the 
letters A and B. 

2.19 A Digression - Combination tables for sets of transformations. 
Consider the parallelogram shown, Fig. 2.13. The transforma¬ 
tions which can be made which leave it occupying the same space 
in the plane are 

(i) a rotation (R) through 180° about O; 

(ii) no movement (the identity transformation). 

These are known as the symmetries of the parallelogram. 




'Hr 

-N* / 

/ 0 

J 




Fig. 2.13 


Let us consider the four possible combinations of R and I, 
writing IR to mean R followed by I (p. 22 explains why this order 
is necessary). We display these in a table 



I R 

Transformation I 


written first 


(performed R 


second) 



Transformation 
written second 
(performed first) 



26 


MATRICES AND THEIR APPLICATIONS 


in which the entry IR appears in the position marked *, that is in 
the I row and the R column. 

Find the four entries in the table, in terms of the single transfor¬ 
mation which would achieve the same final result as the two 
transformations performed successively. For example, IR = R, 
so the entry * is R. 

Make a similar table for the addition of the integers 0, 1 modulo 
2, and note that it has the same pattern as this one. 

2.20 Give a list of the four symmetries of the rectangle, using the 
letters of Fig. 2.14: write the matrix corresponding to each 
symmetry. Form the table of combinations of these symmetries. 



If the reader is familiar with finite arithmetics, find among them 
a combination table which has the same pattern as this one. 
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3.1 The inverse matrix 


The matrix which gives an enlargement of a figure by a factor of k is 

(o *) =Asay - 


Is there a matrix which will restore the figure to its original state? 
Such a matrix will be ^ ^ =B which reduces a figure by 1 jk 

times. The effect of A followed by B is to leave the figure unchanged. 


We may also verify that the product of these matrices is I 



A matrix M is called the inverse of the matrix A, if AM = MA = I. 
Thus B is the inverse of A. 


3.2 Inverses of some known matrices 

What is the inverse of the matrix 

The matrix A = ^ ^ ^ reflects points in the y axis. 

The matrix which will restore points to their original position will 
need to reflect in the y axis again; thus the inverse of A is A itself and 



We can show that if M is a left-inverse of A, i.e. MA = I, then M 
is also a right inverse, i.e. AM = I. Let Y be a right inverse so AY = I. 

In AY = I premultiply each side by M ; then MAY = MI and hence 

MAY = M. 

But since MA = I, and MAY «= (MA)Y by the associative law, IY = M 
thus Y = M and the left and right inverses are the same. 

27 
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The matrix and its inverse thus commute. This can be seen geo¬ 
metrically for if A followed by M restores points to their original posi¬ 
tion, so M followed by A will keep points where they are; the two 
operations are inverse. M can be written as A -1 , i.e. A -1 A = I = AA -1 . 


? 


3.3 Inverse of the general 2x2 matrix 

What is the inverse of the matrix A = ^ ^ 

Here we cannot deduce the inverse from geometrical considerations. 

Let the inverse M of A be ( ^ M . 

\ r V 


Thenas (r ?)(s 2 ) = (o l) gives 


3 p -+- 5q = 1 
p + 2q = 0 
3r+ 5^ = 0 
r + 2s = \ 

by solving the equations. 


p=2 


we obtain ^ ! 

r = -5 

5 = 3 


M is thus 


(-5 3 ) 

[check that (J “*) (5 2 ) = (o l) 
and also (* *) ( J °)] • 


Exercise 

3.1 Find the inverse matrix of 

(i) 


(4 , ); 00 ( 3 _>)> (iO (3 >). 

Is any pattern becoming noticeable? 

To find the inverse of the general 2x2 matrix A = ^ Let the 
inverse of A be M = ^ ^ and it is required to find p , q> r , s in terms 

Of ay by Cy d. 

Fr ° m (r:)(: d)=(o ?) we ° btain 
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ap+qc = 1 
bp+ qd = 0 
ra + sc =0 
rb+sd = 1. 

By solving these equations p = dj{ad-bc) i q= -b/(ad-bc ), 
r= -cl(ad-bc) and s = a/(ad-be). The inverse matrix A is thus 

1 /{ad - be) ^ ^ . Note the pattern of the elements. 

The expression {ad - be) is known as the determinant of the matrix. 

A is and its determinant is denoted by det A or by det 

Its value {ad - be) is found by multiplying the entries in the principal or 
leading diagonal and subtracting the product of the entries in the other 
diagonal. Note that in the initial examples the determinant value was 1. 


Exercise 

3.2 Write down the values of the determinants of the following 
matrices, then write down their inverses and check your answers 

(i) (2 1\ (ii) /I 2\. (iii) /I 2\ (iv) (3 2\. 

\3 lj’ V3 4j’ \4 3 )• \1 5)’ 

(v) ( C ? S a S ^ n =R. Show the inverse of ItfR -1 ) is such that 
' 7 \ — sin a cos ocj v 7 

the rows of R are the columns of R -1 . R said to be transposed. To 

transpose a matrix A, rows of A become columns of the transposed 

matrix, denoted by A'. Then R _1 = R'. 


3.4 The uniqueness of the inverse 


The matrix inverse 
matrix which is inverse 



3 

-4 



Is there another 


Let A be the matrix and suppose X and Y are inverses of A, then 
XA=AX = Iand YA = AY =1. 


From AX = I 

Y(AX) = Y(I) premultiplying each side by Y 
(YA)X=Y associative law for matrix multiplication 
IX=Y since YA = I 

X=Y. 


The inverse of a matrix is unique. 
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3.5 Matrices with no inverse 


What is the inverse of ^ 4 ^ • As we follow the process for finding 

an inverse as in the previous section, the inverse of the matrix ^ 4 ^ 
will have a multiplying factor of the reciprocal of the determinant value 
of the matrix, ^ 4 ^- This determinant value is 1.4 - 2.2 = 0 and so as 

the reciprocal of 0 is not determined, the matrix has no inverse. 


All matrices of the form 


{ka kb) Wi 


will have zero determinants and thus 


no inverse. The important geometrical significance of a matrix which 
has or does not have an inverse is discussed in the next section. A 
matrix which has no inverse is said to be singular , whilst the matrix 
which has an inverse (i.e. its determinant is not zero) is said to be non¬ 
singular. 


Exercise 

3.3 Show diagramatically the figure into which the unit square is 
transformed by the matrix A = ^| Similarly show the trans¬ 
formation of the unit square by the matrix B = (2 4 ) • What is the 

geometrical difference between the transformations of A and B? 

In the case of B, show all points of the plane are mapped onto 
the line y = 2x. The point (1, 1) goes to (3, 6 ). What other points 
are also mapped onto (3, 6 )? 


3.6 Relations between determinants of matrices 
Exercise 

3.4 To investigate whether the determinant of two matrices A and B 
are related in any way. 

Choose any two 2x2 matrices A and B. 

(i) Find detA and detB; then find AB and detAB. Show 
det AB = det A det B. 

(ii) Similarly find A + B and det (A + B). 

Show that in general det (A + B) =£ det A + det B. 

(iii) Compare det A with det 3A and show det 3 A = 3 2 det A. 
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3.7 The translation matrix 

We have so far concerned ourselves with such basic geometric trans¬ 
formations of the plane as reflection, rotation, enlargement and shear, 
and the combination of these. A notable transformation omitted in our 
list is translation. All 2x2 matrices preserve the origin, but in a trans¬ 
lation the origin will be shifted. 



Fig. 3.1 


The plane is to be translated so that the x coordinates are increased 
by a constant h , and the y coordinates by k . The general point P(x , y) 
in the plane is mapped onto P # (x # , y *) where 

x*=x + h m 

y*=y + k * * * • v h 

We have now to extract the translation matrix from these equations - 
which are linear but have a constant at the end; a situation we have not 
met before. 


Equation (1) can be written in matrix form as 



• ( 2 ) 


giving as our translation matrix T = 



This is a matrix 2x3 


matrix which will be awkward to fit and interpret with our 2x2 
matrices. Matrices which are square and of the same order are easier to 
handle. We cannot reduce our T matrix to a 2 x 2, but can increase it to 
a 3 x 3 as 


1 0 h 
0 1 k 
0 0 1 


/x*\ A 0 h\/x\ 

The form (2) then has to become I y* )=( 0 1 & 11 y) and the 

Vl / Vo 0 1/Vi/ 
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equations are 

x* =x + h 
y* =y + k 

1 = 1 

where the last equation is there to see things right! 

Since T is a 3 x 3 matrix, its combination with 2x2 matrices can 
only be achieved if we increase the order of our 2x2 to 3x3. We do 
this by use of the 1 in the bottom right hand corner: e.g. the rotation 
matrix 

( x /cos a - sin a 0\ 

cos a - sin a\ , / . n \ 

becomes ( sin a cos a U I . 
sin a cos a/ \ q 0 U 

The matrix which translates the axes so that points of the plane 

. /! 0 "»\ 

remain fixed and the new origin is (r, s) is I 0 1 - s ). 

VO 0 1 / 


Exercises 

3.5 Find the matrix which will rotate the plane about the origin 
through an angle a, followed by a translation of the plane through 
values (h, k). Give some value to a, h and k and check geometrically 
that the result is correct. 

3.6 Ta is a translation of 2 units in the x axis direction and is a 
translation of 3 units in the y axis direction. and are reflec¬ 
tions on the x and y axes respectively; and R is a rotation through 
90° anticlockwise about the origin. P is the point (a, b). Find the 
coordinates of P 19 P 2 , P& P 4 formed from the point P under the 
following transformations 

(i) P x from P by T*. followed by M r 

(ii) P 2 from P by followed by M r 

(iii) P 3 from P by R followed by T^. 

(iv) P 4 from P by R followed by M v . 

3.7 Find the matrix which will reflect the plane in the line y = mx + c. 
(Hint: translate the axes to (0, c) followed by the matrix to reflect 
the plane in the line y = mx , followed by translation of the axes to 
(0, -c).) 
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3.8 3x3 matrices as transformations of space 

As yet we have considered 2x2 matrices and in the main we shall 
concern ourselves with matrices of this small size. The only reason for 
doing this is simplicity! It is easy to work with 2x2 matrices and the 
methods and ideas used can, in general, be extended to higher orders of 
matrices. We have with 2x2 matrices, transformation of the plane and 
it is straightforward to extend these methods to 3 x 3 matrices as 
transformation of space: we may even extend these to 4 x 4 matrices 
as transformation in 4 dimensions. 


y 


Fig. 3.2 

In 3 dimensions we have the three mutually perpendicular axes 
x, y and z. If we stretch space by a factor of a in the x axis direction, a 
factor of b in the jy axis direction and a factor of c in the z axis direction, 
then a point P(x, y , z) is mapped onto the point P*(x* y y* y z*) where 

x* =ax 
y* — by 
z*=cz. 

The stretch matrix is thus 

a 0 0\ 

0 b 01. 

0 0 J 

Our cube in Fig. 3.2 becomes a cuboid. 

Exercises 

3.8 Find the 3x3 matrices which transform space such that 

(i) reflection takes place in the x-y axes plane 

(ii) enlargement by a factor of k 
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(iii) a shear is performed with a factor of a in the x axis direction 
and a factor of b in the y axis direction (a diagram to see 
the result of the transformation will help before proceeding 
to the equations). 

3.9 Describe geometrically the transformation of space by the 
following matrices: 

(i) /I 0 0\ (ii) /a 0 0\ (iii) /I 0 0\ 

0 -1 0 0 10 ) 01 0 

Vo oi/ Vo o v Vo o v 

(iv) /I 0 0 

0 1 0 

Vi i i 


3.9 Inverses of matrices of higher order than 2x2 

Again we have only considered inverses of 2 x 2 matrices. It is possible 
to find the inverse of a 3 x 3 matrix by solving sets of equations using 
the same idea as used in Page 28. It is instructive to do this (although 
the equations are heavier) and the 3x3 determinant of a matrix is 
obtained. However, an easier method of finding inverses of higher order 
matrices will be given in the next chapter. 
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4.1 Simultaneous linear equations 

The solution of two simultaneous linear equations 

3x + 2y = 8 
4x + 3y = 11 

is as follows. 


These equations can be written as ^ ^ = (l l) ’ t ^ iat ls ^ ^ ls 

ix ^ ^ and k is the vector 


the matrix 


then the set of linear 


equations can be written as Ax = k. 

Premultiply each side by A -1 ; then 

A -1 Ax = A“ 1 k, 
Ix = A _1 k, 
x-A^k. 


To solve the given set of equations we need to premultiply ^ ^ by the 

G?)- 

is ( - 4 s)’ 80 


inverse of 

The inverse of A 


G)-(4 D(n)' 


Thus 




or 


x = 2 

y = \' 


There are, of course, more direct methods of solving these equations 
We give the matrix method here as an example of a method which 
generalises immediately to any number of equations. 
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Exercises 

4.1 By matrix methods solve these sets of equations, and check your 
answers by substitution: 

(a) 4# + 5y = 23 (b)3x-y=U (c) 2x + 3y = 4 

x + 2y = $ 2x + 5y= -4 - 6# + 3y = 8 


4.2 By matrix methods solve the sets of equations given by Ax = h 
where A = ^ ^ ^ and h takes in succession the values 

(,, ($ (ii) (Is) ( 13 ) W (-!) 

(vi) own choice of values. 


4.2 The nature of solutions 

As we have shown, the solution to the set of equations 

3x + 2y = 8 
4x + 3y = 1 

is# = 2, andjy = 1. 

What would be the solution if the equations were homogeneous; that is 


3x + 2y = 0 
4x + 3y = 0? 

The only solution is x=y = 0. 

The matrix ^ singular (since the determinant is zero) so this 

matrix has no inverse. What is the solution to the set of equations 

2x + 4y = 6 
x + 2y = 3? 


The solutions are x = 3 - 2t, y = t for any real number t ; an infinite set of 
solutions. 

What is the solution if these equations are homogeneous, i.e. h = 0? 
This time we have solutions x = 2 £, y = - 1, again an infinite set of 
solutions. 


The equations 


have no solution, 


2x + 4y = 6 
x + 2y = 2. 
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The nature of the solutions of a set of linear equations can be sum¬ 
marised as 



non-homogeneous 

(h#0) 

homogeneous 

(h=0) 

(i) non singular 

unique 

solutions = 0 

matrix 
(det A t* 0) 

solution 


(ii) singular matrix 

many solutions or no 

many solutions 

(det A=0) 

solutions (depending 
on h values) 



Exercise 

4.3 Make sets of equations in each of the above five cases and show 
the graphical interpretation for solutions. 

The part of this work which will be required later is that in a set of 
homogeneous equations, if solutions other than zero are required the 
equations must be such that det A = 0, i.e. the matrix must be singular. 


4.3 Inverse of a 3 x 3 matrix 

The set of simultaneous linear equations 

#-3;+ 2# = 6 'j 

3x-2y+z = 7 > . . . . (1) 

x+ 3y-z =-4 j 

is solved conventionally by multiplying the first equation by - 3 and 
adding the result to the second equation; multiplying the first equation 
by - 1 and adding to the last equation to give 

x -y + 2z = 6 ■) 

y-5z=-ll\ ... .(2) 

4;y-3*= -10 J 

The set of equations (2) is equivalent to the set (1). 

We continue to solve the equations by multiplying the second equa¬ 
tion by - 4 and adding the result to the last equation, to give 

x-y -\-2z-6 'j 

>>-5#=-lli . . . .(3) 

17#=34 J 
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From here we obtain the value of z and, by substitution, the values 
first y and then x as 


x = l 

y =-1 

z = 2 



. (4). 


of 


This method of solution, reducing the set of equations to a triangular 
form and ultimately to a diagonal form (4) gives an indication of the 
process to be used with matrices. (It also provides a systematic method 
for solving sets of linear equations.) The set of equations (1) in matrix 
form is 

(1 : i j)0=(J) « - 


and the solution set (equations (4)) is 


/I 0 O 
0 10 
VO 0 h 




The process of solution can be regarded as the method by which 

A -1 2 \ A 0 0 \ 

I 3 - 2 1 J is changed into 10 1 0 1 

Vl 3 -V Vo 0 1/ 

by a process called row reduction - a sequence of operations performed 
on the rows of the matrix. 

/ 1 0 °\ 

Consider the matrix E^l -3 1 0 1 with the matrix 

V-l 0 1/ 

A -1 2 \ 

A = ( 3 -2 1 ), noting that is basically the 3x3 identity matrix 

Vl 3-1/ 

( \ 0 0 \ 

I 0 1 0 J with the second two elements of the first column of A used 

Vo 0 1/ 

with changed sign. 

/ 1 0 Owl -1 2\ A -1 2\ 

ThenE 1 A=( -3 1 0 3 -2 1= 0 1 -5 givesthe 

V-l 0 l/ Vl 3 -1/ Vo 4-3/ 

left-hand set of equations (2). The matrix E x has the property of per¬ 
forming the row reduction required in solving the equations since the 
row (- 3 1 0) in E 2 , when matrix multiplied with A multiplies row 1 
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of A by - 3 and adds it to row 2 of A - the method used at the beginning. 
By using E x with h as E x h the right-hand side of equations (2) is obtained. 


Continuing to solve we would choose E 2 as I 0 

Vo 

f l -l 2\ 

before E X A = 10 1 - 5 1 from above; so 

Vo 4-3/ 

/I 0 Owl -1 2\ /I 

E 2 E 1 A = (0 1 0(0 1 -5 = 0 

Vo -4 1 /Vo 4-3/ Vo 


0 0 \ 

1 0 ) and use it 

-4 1/ 


-1 

1 

0 



Thus E 2 E 1 Ax = E 2 E 1 h and gives the set of equations (3). At this stage 
we solved the equations by finding z and substituting. We could 
continue using matrices. 



/ 

1 0 

On 




Choose Eg 

as (i 

0 1 

0 

) then 



V 

0 0 

iV 

1 



A 0 

°\ 

A 

-l 

2 \ , 

/I 

-1 

(0 1 

0 ) 

(° 

l 

-5) = 

0 

1 

Vo o 

& 

Vo 

0 

17/ 

Vo 

0 


We now chose E A as 


the 5). Then 


-5 ) = (EgEgEjA). 
V 


1 0 0 \ 

0 1 5 j (note the basic identity matrix with 

0 0 1 / 


/I 0 Owl -1 2\ /I 

0 1 5 (0 1 -5=0 

Vo 0 1 / Vo 0 1 / Vo 


-1 

1 

0 


=E 4 E 3 E 2 E 1 A 


and the right-hand side would be E^gE^jh and would give the value 
of y, 

A 

Our final choice for E 5 to obtain x would be ( 0 


Vo 


-2 

0 

1 


. Then 


z 1 1 

~ 2 \ 

/i 

-1 2 \ 

A 0 

°\ 

(° 1 

°) 

(° 

1 o) = 

=(o 1 

0 ) 

Vo o 

1 / 

Vo 

0 1 / 

Vo o 

1 / 


= E 5 E 4 E 3 E 2 E 1 A. 


Thus by a series of matrices EgE^E^ we have changed the original 


A 

-1 2 

\ f l 

0 

°\ 

matrix ( 3 

-2 1 

) into [ 0 

1 

°) 

Vi 

3 -1- 

/ Vo 

0 

i/ 


inverse of A (as E 5 E 4 E 3 E 2 E 1 A = I). 
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1 1 -2 

e^ejea^o 1 0 

0 0 1 


1 

0 

°\ 

/I 

0 

°\ 

A 

0 

O' 

0 

1 

5) 

(° 

1 

0 ) 

(° 

1 

0 

^0 

0 

\) 

Vo 

0 

& 

Vo 

-4 

b 


1 0 0 
-3 1 
-1 0 


?) 


“17 

o 

17 

17 \ 

4 

17 

3 

“17 

Tf =Esay. 

11 

4 

JL/ 

17 

“17“ 



Check that EA=I and also Eh gives the solution for x 9 y and z. 

At present this simply appears as a viable but somewhat tedious 
method of finding an inverse - but now understood as a process it can be 
presented differently and ‘speeded-up’. 

1-1 2 \ 

3 - 2 1 ]. We did this by 

13-1/ 

which eventually changed the 


Our task is to find the inverse of A = 


using matrices E 1? E a , E 3 , E 4 and E 5 


1 


0 


0 

matrix to ( 0 1 0 }. The effect of the E matrices is to perform row 

^0 0 V 

reductions. We can compute the effect of the E matrices as we proceed 
so that at the end we shall not have to multiply them out to obtain the 
inverse. To be able to start the computation we use the matrix 

f\ 0 0\ 


r u u \ 
=(0 1 °). 
Vo o 1 / 



A 


i 


(1) /I 

-1 

2 

1 0 

0 

(3 

-2 

1 

0 1 

0 

Vi 

3 

-1 

0 0 

1- 



/ 1 

0 0\ 

/l 

Now 

EjA = 

1 ~ 3 

1 0 1 

3 



V_i 

0 1/ 

Vi 

/ i 

0 C 

)\ 



E*I= -3 

1 c 

) ) . We may put i 

V-l 

0 1 

./ 



and I as 





(2) /I 

-1 

2 

1 

0 

0 

1 

-5 

-3 

1 

Vo 

4 

-3 

-1 

0 


-i 

-2 

3 


?)-(i 


-l 

l 

4 


2 

-5 

-3 


and 


°\ 

;)■ 
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The effect of using E x on A and I is to perforin row reductions, and 
these reductions are performed just as we did originally on the equa¬ 
tions. In (1) row 1 times - 3 to add to row 2, row 1 multiplied by - 1 to 
add to row 3 to give rows 2 and 3 in (2). 


A 0 Ox 

Proceeding to use E 2 =l 0 1 0 1 with EjA and E]I we obtain 

VO -4 V 


A -l 2\ 
E 2 E 1 A= 0 1 -5) 

Vo 0 17/ 


results together as 

(3)/l -1 2 

0 1-5 

Vo 0 17 


and E 2 E 1 I = (^ -3 

V 11 


1 0 0 \ 
-3 1 0 

11 -4 1/ 


0 0 \ 

1 0 ). Putting these 
-4 1/ 


we see row reduction has again been used. The effect of E 2 is to take 
row 2 in (2), multiply it by - 4 and add it to row 3 to give row 3 in (3). 
By doing this reduction to the whole row in (2) we are computing 
E 2 EjI on the right-hand side and hence the E matrices giving the inverse 
of A are being computed as we proceed. 

E 3 E 2 E X A and EgE^I are computed by dividing row 3 in (3) through¬ 
out by 17, to give row 3 in (4). 


(4) 

e 4 


1 -1 2 
0 1 -5 

0 0 1 


is 


1 

-3 

li 

17 


0 

1 

4 
' 17 


1 0 0 \ 

0 1 5 ), which has the effect of taking row 3 of (4) multi- 

0 0 V 


plying it by 5 and adding it to row 2 to give row 2 in (5). Thus for (5) 
we have E 4 E 3 E 2 E 1 A and E 4 E 3 E 2 E X I as 


1/1 -1 

2 1 

1 

0 

°\ 

f° 1 

0 1 

4 

17 

3 

“17 

fV). 

Vo 0 

1 

11 

17 

4 

“17 


/I 

1 

" 2 \ 



for E 5 = ( 0 

1 

0 ) 

which 

adds row 1 to row 2 and multiplies 

Vo 

0 

1 / 




row 3 by - 2 and adds it to row 1 to give rows 1 and 2 in (6), giving 
E 5 E 4 E 3 E 2 E 1 A and E 5 E 4 E 3 E 2 E X I computed as 


(6) i 

/i 

0 

0 

1 

“17" 

5 

17 

1 3? \ 

0 

1 

0 

4 

17 

3 

“17 

tt) 


Vo 

0 

1 

11 

17 

4 

“17 

tV 
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We have the inverse of A on the right-hand side. The clue to the process 
is the changing of A by row reduction into the identity matrix, perform¬ 
ing the same reductions on the right-hand side. Obtaining the identity 
matrix on the left-hand side is a systematic process which is as follows: 

( a b c\ 
d e /) then 
g h i) 

(i) obtain a 1 in place a - divide top row throughout by a 

(ii) obtain zeros on places d and g - multiply row 1 appropriately 
and add 

(iii) obtain a 1 in place e 

(iv) obtain zeros in places b and h 

(v) obtain 1 in place i 

(vi) obtain zeros in places c and/. 

(This reduces some of the work we have done above - more could be 
done in (3) and (5) to obtain zeros.) 


Example 

( 2 1 

To find the inverse of A = 1 4 -1 

V-l 1 


A 


I 


( 2 1 
V-l 1 

/ 1 i 

4 -1 
V-l l 

( 1 * 
0 -3 

V 0 | 

i 4 
o 1 
o i 

( 1 0 
0 1 

Vo o 


-3 

1 

0 

°\ 

-2 

0 

1 

°) 

2 

0 

0 

1 / 

3 ! 

~ 2 

1 

2 

0 

°\ 

-2 

0 

1 

0 

2 

0 

0 

1 / 

3 

~ 2 

1 

2 

0 

°\ 

4 

-2 

1 

0 

1 

2 

1 

2 

0 

\) 

3 

“ 2 

1 

2 

0 

°\ 

4 
— 3 

2 

3 

1 

~ 3 

°) 

1 

2 

1 

2 

0 


5 

6 

1 

6 

1 

6 

°\ 

4 
~ 3 

2 

3 

_JL 

3 

°) 

5 

2 

1 

“ 2 

1 

2 

1 / 



Instructions for working from 
one set to the next 

Divide row 1 by 2 

row 1 x - 4 add to row 2 
row lxl, add to row 3 

divide row 2 by -3 

row 2 x - \ add to row 1 
row 2 x - f add to row 3 

divide row 3 by f 
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( 

( 


1 

0 

0 

1 

0 

0 


0 -t 
1 -I 
0 1 

0 0 
1 0 
0 1 


JL 

6 

2 

3 

JL 

5 

0 

2 

5 

1 . 

5 


1 

6 

JL 

3 

1_ 

5 


O' 

0 



row 3 x f add to row 2 
row 3 x f add to row 1 


This is still an arithmetical burden, but it presents a systematic process 
(which can be programmed for a computer) of finding the inverse of 
any (square) matrix. Another method involving determinants is available 
for 3 x 3 matrices, but it is perhaps best to keep to the given standard 
method for all cases. 


Exercises 

4.4 Find the inverse of the matrices 




0 

1 

3 



4.5 Choose a 2 x 2 matrix and find its inverse by choosing the suitable 
E matrices. Now look at the geometrical aspect. Use your matrix to 
transform the unit square. With each of the E matrices show diagrama- 
tically how they systematically restore the transformed figure back to the 
unit square. 
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THE EIGENVECTOR 
5.1 Invariants in a mapping 

In this section we shall consider what properties of a figure remain 
unchanged when a transformation is made. The matrix to rotate the 

plane through an angle a about O is ^ ^ i.e. if y*) 

is the new point and P(x , y) the original point then 



Fig. 5.1 

Thus for the unit square in the diagram, Fig. 5.1 

A( 1, 0)->^*(cos a, sin a) 

B( 1, l)-^P # (cos a - sin a, cos a + sin a) 

C(0,1)->C*( - sin a, cos a). 

Although the points A> B, C have moved in the transformation the 
figure formed by OABC is still a square; the lengths of OA , BA etc. and 
the angles between the lines remain unaltered. So the expressions for 
these entities in terms of the new coordinates must be expressions which 
yield the same value as they did in the old coordinates. For example, 
OA*= s /(cos 2 a + sin 2 a) = l =OA. These are called invariants of the 
mapping. 
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Exercise 


5.1 Find the length B*C * and the angle, lOA*B* in the new co¬ 
ordinates and show that their values are the same as in the old 
coordinates. 

On the other hand, the direction of the line OA is represented by the 
vector (1,0) in the old position and by the vector (cos a, sin a) in the 
new position; so direction is not an invariant. 

Does the matrix ( C .° S a Sln ) provide an y other invariants in the 

\sm a cos ocj J 

geometry of its mapping? What of area? The area of the square remains 
the same, so here is another invariant. 


If we consider the general matrix 

of its mapping? 

Again consider the unit square OABC. 



can we find the invariants 



Under this transformation (Fig. 5.2) 

o->o 

A( 1, 0 )->A*(a, c) 

B{\,\)^B*(a + b,c + d) 

C(0, 1 )->C*(6, d). 

Lengths and angles are now not invariant. O remains at O. What of 
area? - the area of OA*B*C* is not the same as that of OABC , so area 
is not an invariant. 


Exercises 


5.2 


(a) Using the matrix 



on the unit square OABC , show that 


the transformed figure OA*B*C * is a parallelogram. 
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(b) Show the area of the parallelogram OA*JB*C* is (ad-be) 
(the determinant of the matrix). (Express the area as the sum or 
difference of triangles and trapezia as in Fig. 5.3.) 



(c) Deduce that the condition for the area to be an invariant of 
the mapping is (ad-be) = ± 1. 

5.3 (a) For the following matrices, find the value of their determinants 

(i) (2 1\ (ii) (3 2\ (iii) (2 4\ (iv) (l 3\ 

\5 1) ^5 4 ) \3 6) \3 6) 

(v) / cos 2 a - sin 2a\ 

\sin 2a cos 2a J 

(b) Which will have area as an invariant of the mapping? 

(c) Show geometrically the effect on the unit square of the 
mapping represented by these matrices, explaining the relationship 
between the determinant of the matrix and the change in area 
produced by the mapping. Note that thq determinant gives the 
numerical value of the area. A minus sign implies a reverse order 
of the vertices, i.e. a reflection in the mapping. 


5.2 The eigenvector as an invariant of a matrix mapping 

Besides lengths, angles and area, is there any other invariant which a 
matrix mapping may have? The following investigation leads to an 
important invariant possessed by all mappings produced by matrices. 


vector 


Consider the effect of using the matrix A = ^ j 5^ to mu ^Pty the 
(l) * We ^ ave 

(-S sXKMO 

so that the result is to multiply by the scalar 2. 
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This is an interesting and possibly significant situation, particularly 
when we consider the representation of the vectors as displacements, 
from the origin, in a plane, and the matrix as a transformation of the 
vectors. The situation is then that each of these vectors (u and v in 
Fig. 5.4) is transformed into a vector in the same direction as itself; 
u->2u and v->4v. 

The question now arises whether the matrix has the same effect on 

any other vectors, or whether and are special. Try w for 

example; this and its transform are shown in the diagram; they are not 
in the same direction. 

Let us attempt systematically to find whether there are any other 


vectors ( J which behave with this matrix like 

\yj 


DGM;)’ 


. We must 


where A is a scalar 


-h 

-x + 5y=Xyj 

(l-A)* + 3;y=0\ m 

-x + (5-A)y = 0J ■ ■ ■ ■ V>- 

From the consideration on page 37, if there are to be solutions to these 
equations other than x =y = 0 we must have 


0-1 !-»)■ 
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Hence 

(1 -A)(5 - A) + 3 = 0, 

A 2 - 6A + 8 = 0, 

(A - 2)(A - 4) = 0 

A = 2 or 4. 

(These are the two multiplying factors we found previously.) 

To find whether any other vectors exist we put each of these 

values of A in turn, back into the equations (1). Do this, first with A = 2 
in both equations. What do you notice about the two equations? And 

what values of are solutions? Are there any others besides ^j^? 
What about or Now put A = 4 into both equations and again 


find the solutions for . State your conclusions about the existence of 

further vectors which behave like and ^ . 

We now give these special vectors a name. 

If A is any matrix, x a vector and A a scalar, such that Ax = Ax, then x 
is called an eigenvector of A, with eigenvalue A. 

Example 

Find the eigenvalues and eigenvectors of the matrix Q ^ • 


Let A = ^ 2 J anc ^ ^ et x an eigenvector of A with eigenvalue A, 

then Ax = Ax and ( 1 2\ ( x\ . ( x\ 

[l -2){ y)-\y)’ 
x + 2y =\x 
2x-2y =A y, 

and (1-A)^ + 2 j = 0 > 1 m 

2x-h(-2-\)y = Of ' ' * * (L) ‘ 

For equations (1) to have solution other than x = 0 =y then 

“"G"* -2-a)-° 

(1 -A)(-2-A)-4 = 0, 

A 2 + A-6=0, 

(A-2)(A + 3) = 0 

therefore A = 2 or - 3, the eigenvalues. 
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When A = 2, the equation (1) becomes -# + 2y = 0 and 2#-4y = 0, 
both equations being the same; x-2y = 0. This is the eigenvector line 

corresponding to the eigenvalue of 2. Particular eigenvectors are 
found by finding values of x and y which satisfy x - 2y = 0; e.g. ; 
^ 2 ^ ’> et c. Note that for any such vector, say > the matrix A has the 


effect of multiplying by 

fS 




thus 


(1 2 

\2 -2 


= " =2 


When A = -3, the equations (1) become 4-x + 2y = 0 and 2x+y = 0. 
Hence the eigenvector line is 2x + y = 0, with an eigenvector ^ ^ • 


5.3 Finding the eigenvectors of a matrix 

We must summarise the above work by dealing with the general case. 

The eigenvectors of the matrix A = ^ are found by considering 

the mapping of P(x y y) onto P # (A,v, A y ); that is ^ ^ = ^(^j • The 

set of linear equations giving this mapping are 

(a-X)x + by = Cfi m 

cx + (d-X)y = 0J ’ ‘ ‘ 

or A x x — 0, where A x = 


(a-X b \ 
\c d-XJ' 


The condition that equations (1) have solutions other than x = 0 =y is 
det A x = 0; that is 

*(r A i. 

or (a-\)(d-X)-bc = 0 .... (2). 


?-a)~ 


0 


Equation (2) is known as the characteristic equation of the matrix. 
The values of A from (2) are substituted into (1) to give the eigenvector 
line - an equation of a line through the origin such that any point on the 
line is mapped (by the matrix) onto another point on the same line. 

Since 

IH IK. 5K* 

the characteristic equation can be written as det (A - AI) =0 
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Exercises 

5.4 Find the eigenvalues, eigenvector lines and eigenvectors of the 
matrices 



5.5 Find the eigenvalues of ^ and 

5.6 Write and simplify the characteristic equation of the matrix 



Deduce what you can about the sum and product of the eigen¬ 
values. 

5.7 Is it possible for a real matrix to have complex eigenvalues? Try 
to make one up. 

5.8 Is it possible for a real matrix to have a single repeated eigenvalue? 
Try to make up some which have. 

5.9 What can you say about the eigenvalues of a diagonal matrix? 
(One with zeros everywhere except on the principal diagonal; 
write down one and try it.) 

5.10 What are the eigenvalues of a matrix of the form kl ? 

5.11 What are the eigenvalues of the matrix 

5.12 Find and note the eigenvalues of the following matrices: 

G9G9’ (J-9G-9- 

5.13 Find the eigenvector lines of the matrices 

G 9- (o -!)• (-? -9- (-? -9- 

5.14 Investigate such further matrices as you need to, to draw con¬ 
clusions about the eigenvectors of (a) matrices with real, distinct 
eigenvalues; (b) matrices with complex eigenvalues; (c) matrices of 
the form kl; (d) matrices with a single repeated eigenvalue; 
(e) diagonal matrices. 
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5.4 The eigenvectors of a 3 x 3 matrix 

(2 " 2 3 \ 

Let A = I 1 1 1 1, then the eigenvector of A, x is given when 

Vl 3-1/ 

Ax = Ax, where A is the corresponding eigenvalue. With A as a 2x2 
matrix, the eigenvector is two-dimensional; now A is a 3 x 3 matrix and 
x is a three dimensional vector. 

Hence 



2x -2y + 3z = Ax 
x +y + z = Ay 
x + 3y -z=Az, 

(2 - A)# - 2y + 3z = (h 
x + (l-X)y + z = 0\ . . . .(1). 

x + 3y + {- \ -X)z = 0 j 

As before, equations (1) have solutions other than x = 0 =y = z if 

/2 - A -2 3 \ 

det 1 1 -A 1 =0. 

Vl 3 -1-A/ 

The cubic equation for A, the characteristic equation for A, may be 
obtained by expanding the third order determinant, or we may obtain it 
by eliminating x, from equations (1) and then using a second order 
determinant. This latter process is as follows 

(2-A)x-2y+ 3z = 0 . . . . (a) 

# + (l-A)y+ £ = 0 . . . .(b) 

x + 3y + (-l- A)# = 0 . . . .(c) 

(2 - A) (b)-(a) gives 

(l-A)(2-A)^ + 2y + (2-A)«-3« = 0 

that is 

(A 2 -3A + % + (-l-A> = 0 . . . .(d) 

(b) - (c) gives 

(-2-A)y + (2+A)* = 0 . . . .(e) 
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From equations (d) and (e) we have solutions other than y = 0 = z 

whendet( A2 :f_ + A 4 _ 2 + a) = ° 

(A 2 - 3A + 4)(2 + A) - (- 2 - A)( -1 - A) = 0, 

A 3 - 2A 2 - 5A + 6 = 0, 

(A - 1)(A + 2)(A - 3) = 0. 


The eigenvalues of A are 1, - 2, and 3. 
When A = 1, in equations (1) 


x - 2y + 3z = 0 . 

• • • W. 

x +#=0 . 

. . . (ii), 

+ 

i 

to 

II 

o 

. . . (iii). 

From (ii) x= -z, so (i) becomes 

-2y + 2z = 0 and (iii) becomes 

3y - 3z = 0. 


Both of these equations givejy = z. 


The eigenvector line is given by 

xj - l=yjl =#/l, which is a line 


through the origin. A particular eigenvector is ^ 1^ . Note that 



When A= -2, the eigenvector line is x/11 =y /1 = z/ - 14, an eigen¬ 


vector is 


n \ 

1 1; and when A = 3 the eigenvector line is x/1 =y/ 1=2/1. 
-14/ 


Exercise 

5.15 Find the eigenvalues and eigenvectors of the matrices 


(0 


/i 

0 

_1 \ 

(ii) / 1 

1 

i 

2 

i 

- 1 

2 

V2 

2 

3/ 

V o 

1 



5.5 Transpose of a matrix 

The transpose A' of a matrix A is the matrix formed from A by writing 
he rows of A as columns; the wth row of A becomes the nth. column 
of A'. 
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Thus the transpose (A') of the matrix A = 




If a matrix is square its transpose is a reflection in the leading diagonal; 


A = 




5.6 Transpose of a product 
Exercise 

5,16 A = (* j) andB = ^ _ 2 ^ find AB, A'B',B'A',(AB)'. 

Is there any relationship between A 7 , B' and (AB)'? 

Using 2x2 general matrices, prove the relationship between 
A', B' and (AB)' is (AB)' =B'A'. 


5.7 Orthogonal vectors 

u x 
U 2 

v = ^ x ^ the product u'v is (u v m 2 ) ^ x ^ = (u x v x + w 2 z; 2 ) - a scalar; in fact 


^ is a column vector, then u' = (u v w 2 ) is a row vector. When 



the ‘inner product* u . v of the vectors u, v. If u = 



and 



then u' . v = 0, and as the gradient of u is 3 and that of v is - u and v 
are perpendicular. If u' . v = 0 the vectors u and v will be perpendicular 
or orthogonal. 


5.8 Symmetric matrices 

If a matrix A is such that A = A', then A is said to be symmetric. A 
symmetric matrix must be square. 

Exercise 

5.17 Find the eigenvalues and show that the eigenvector lines of the 
symmetric matrix ^ ^ are orthogonal. Show that the general 

symmetric matrix M = ^ ^ always has real and distinct eigen¬ 

values and orthogonal eigenvectors. 
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5.9 Inverse of a matrix - another method 


The characteristic equation of the matrix A = 



can be verified 


as A 2 - 8A + 7 = 0. Replace A by A, then 


A»-8A + 7.Q *)‘-8g J)+7(J J) 

_/17 32\_/24 32\ (1 0\_/0 0\ 
\16 33/ \16 40/ + \0 7J~\0 0 )’ 


A satisfies its own characteristic equation. 


The inverse A -1 of A can now be found as follows: A 2 - 8A = - 71, 
A(A - 81) = - 71 (Distributive Law). 

Premultiply each side by A -1 , then A _1 A(A - 81) = A -1 . - 71, that is 
A-81= -7A -1 

hpnrp 

A-~«A-8D--i[g t)-(o 2)] 

=-K1 D- 

The crux of this method is the fact that the matrix A satisfies its own 
characteristic equation. Will this always be true? Write down other 
2x2 matrices (or the general 2x2 matrix) and see if they satisfy their 
characteristic equations; then find their inverses by this method. It is 
possible by this process to find the inverse of 3 x 3 matrices. (The 
general theorem that every square matrix satisfies its characteristic 
equation is proved in Birkhoff and MacLane, A Survey of Modern 
Algebra , 1953 edition, p. 320.) 


Exercise 

5.18 Show that the characteristic equation of the matrix 

A 0 3\ 

M=( 1 1 0 
V2 1 4/ 

is A 3 - 6A 2 + 3A -1 = 0. Replace A by M and establish that 
M 3 - 6M 2 + 3M -1 = 0. 

Multiply through by M _1 , showing M -1 = M 2 - 6M + 31 and hence 
find M -1 . 
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APPLICATIONS OF EIGENVECTORS 
TO LINEAR MAPPINGS 

6.1 Introduction 

This section consists of an extended investigation in which eigen¬ 
vectors play the crucial role. It is not, however, essential for under¬ 
standing the later parts of the book. 

The investigation consists of trying to discover what transformation 
of the plane a given matrix will perform. It may be an enlargement, 
reflection, rotation, shear or projection or combination of these. 

The general basis of the method will be to rotate the axes of co¬ 
ordinates so that they lie along the axes of the transformation (for 
example in the case of a reflection, along and perpendicular to the 
mirror), thus reducing the matrix of the transformation to an easily 
recognisable form. In seeking the axes of the transformation, the eigen¬ 
vector of the matrix will have an important role. 


6.2 Transformation matrices and their eigenvectors 


The reflection matrix 



which reflects points in the x axis, 


leaves points on the x axis where they are and reflects points of the y 
axis onto other points of the y axis. So the axes are eigenvectors of this 
reflection matrix; and since ( x , o)->(#, o), while (0, j>)->(0, - y ), we see 
that the eigenvalues are 1 and - 1 respectively for the two eigenvectors. 

The relation between the different kinds of mapping and their eigen¬ 
vectors is shown in Table 6.1 (p. 56). 


6.3 To find the mapping of a given matrix 

We shall now consider how far this set is sufficient for expressing any 
given mapping. We may note straight away that we do not need all three 
of D, E and L, since by taking k = b and /= ajb , we have D =EL=LE. 

As a first example of what we propose to do, consider the matrix 
cos 2/3 sin 2/3\ 
sin 2/3 — cos 2/3 J ' 
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M 


-c 


Its eigenvalues are readily shown to be +1, 
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Matrix 

description 

deter¬ 

minant 

eigen¬ 

values 

eigenvectors 

(a) D = 

(S » 

double stretch 

ab 

a y b 

x axis, y axis 

(b)F = 

C. -?) 

reflection in x 
axis 

-1 

1, "I 

x axis, y axis 

(c) E = 

(S» 

enlargement 

k 2 

k 

any vector 

(d) L = 

(o!) 

linear stretch in x 
axis direction 

l 

ly 1 

x axis, y axis 

(e) S = 

(11) 

shear in x axis 
direction 

1 

1 

x axis 

(f) P = 

Co!) 

projection parallel 
y axis into x axis 

0 

1,0 

x axis, y axis 

(g) H- 

f cos a - sin a\ 
\sin a cos a ) 

rotation of points 
through a 

1 

cos a 

dz 

not real 





i sin a 



(«, b y ky l are >0) 
Table 6.1 


-1 and its eigenvectors and ^ ^ • (Note that they are 

orthogonal and at an angle ft to the original axes.) This suggests that M 
represents a reflection in the line making an angle /3 with the x axis, 
which can be verified directly but which will also be proved by what 
follows here. Suppose we now consider a rotation of axes through an 
angle j8; the reflection is now in the new x axis, so the matrix M, when 


transformed in the appropriate way, must become ^ which is 

one of the basic set we have chosen. It remains to decide what is the 
appropriate way to transform the matrix of a mapping when the axes 
of coordinates are rotated. Figure 6.1 shows original axes (#, y) y new 
axes (x l9 rotated through an angle j8 from the old, and a typical point 
P which maps into P # . The coordinates of P and P* may be expressed 
with respect to either the old or the new axes. If the matrix of the map¬ 
ping is M with respect to the old axes, we have writing x for (^j and 

Xlf0r (J x*=Mx .... (1). 


The relation between the old and new coordinates of P and P* is 
given by 

x x = Rx 
Xj*=Rx* . 


and 


• ( 2 ) 
. (3) 
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where R = ^ sinjS cos/!) ~ t ^ ie matr ^ x f° r expressing the change in 

the coordinates of points when the axes are rotated through +/?, see 
p. 15. 

Thus we have 

x x * = Rx* from equation (3) 

= RMx from equation (1) 

= RM(R" 1 x 1 ) from equation (2) 

= (RMR _1 )x 1 

which shows that the matrix giving the same mapping as M but with 
respect to the new axes is RMR -1 . 



Since this process is very important in the remainder of this section, 
and is not easy to understand, we give a second explanation of it which 
may help those who find the above difficult. 

We are still referring to Fig. 6.1. 

We need to rotate the axes through p into the x l9 y v axes position, 
keeping the points fixed but referring them to axes rotated through ft. 
If we consider a general point P in this mapping, then under a matrix M, 
P->P* where x* =Mx. The old points in the plane, when the axes are 
rotated through an angle /?, will become a set of points with respect to 
the new axes x Xl y x given by Xj=Rx, where R is the rotation of axes 
matrix. Similarly the transformed points under the matrix M, will have 
coordinates with respect to the x v y x axes given by x x *= Rx*. 

The transformed points taken with respect to the rotated axes are 
given by 

x x * = Rx* 

=R(Mx) (from x* = Mx) 

= RM(R _1 x 1 ) (as x x = Rx so that x = R“ 1 x 1 ) 
=RMR" 1 x 1 . 

So the effect of rotating the axes is that the matrix of the mapping 
changes from M to RMR -1 with respect to the new axes. 

c 
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We now apply this result M->RMR _1 to the matrix we are consider¬ 


ing, 


we have 

RMR -1 = ( C08 J 
\ - sin p 

_ / cos /2 
\ - sin /2 


cos 2/2 sin 2/2\ # 
sin 2/2 - cos 2/S J ’ 


sin /2\ / cos 2/2 sin 2/2\ / cos /2 
cos /2/ \sin 2/2 - cos 2/2/ \sin /2 

sin /2\ 
cos /2/ 


- sin /2\ 
cos /2/ 


/cos 2/2 cos /2 + sin 2/2 sin /2 - cos 2/2 sin /2 + sin 2/2 cos /2\ 
\sin 2/3 cos /2 - cos 2/2 sin /2 - sin 2/3 sin /3 - cos 2/2 cos /2 ) 



sin /2\ /cos /2 sin /2\ 

cos/2/\sin/3 -cos /2/ 

= F as named in the list on p. 56. 


Thus we have proved that when the axes are rotated through an 
angle /2, the mapping produced by the matrix M becomes a reflection 
in the x x axis; so with respect to the old axes it must have been a reflec¬ 
tion in the line making an angle /2 with the x axis. 


If we look at this last relation from another point of view, and trans¬ 
form as follows: 

RMR- 1 = F 

MR“ 1 = R~ 1 F, 

MR 1 R = R 1 FR, 

M = R 1 FR. 


We have achieved our aim of expressing M as a combination of matrices 
chosen from the basic set. 


6.4 Another mapping 

As a second example, consider the mapping produced by the matrix 
K = f j To transform this into a combination of matrices of the 


basic set, we have to find its eigenvectors and then rotate the axes of 
coordinates so that the new axes are in the directions of the eigen¬ 
vectors. The eigenvalues turn out to be 4 and 2, with eigenvectors 


u = 



and v = ^ 


respectively. Figure 6.2 shows the unit square and 
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the parallelogram into which it maps, and the two eigenvectors of the 
mapping. 

We now choose one of the eigenvectors - it does not matter which - 


to become the new x axis (a^). We choose 



This is at 45° with the old 


x axis so we choose axes x l9 y x at 45° to the original x, y axes and use the 
matrix for rotating axes through 45°. 


The matrix R required will be ( C ? S 0 that is 

n \-sm45 cos 45 / 

£(-! {)•-**-£(! -!)■ 


The matrix representing the same mapping with respect to the new 
axes is given by 

““■"-^(-1 !)(-! 5 ) 72(1 "!) 

■ i (-2 *)(! 1 ) 

<tt) 

-(0 “)(o 1) 
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which is a shear of factor 1 followed by a double stretch of factors 4 and 
2; these mappings being taken with respect to the x l9 y x axes. Thus this 
analysis of the matrix has enabled us to determine the nature of the 
mapping. 


If we had taken the eigenvector 



as an axis a similar result 


ensues. 


6.5 Further examples of mappings 

There now follow three worked examples on finding the nature of 
the mapping produced by a matrix. In these, various aspects of using 
the eigenvectors as axes have been developed. The chapter ends with a 
set of miscellaneous examples for the reader. 



The characteristic equation of M is (1 - A)( -1 - A) = 0, from which 
A = ± 1. The equations to give the eigenvectors are 


(1 -X)x -2y — 0 . . . . (1) 

and 

0 . * + ( - 1 -A)ry = 0 . . . .(2). 


When A = 1, from (1) 0 . x - 2y = 0, i.e. x can have any value and y = 0; 
which gives ; i.e. the x axis as an eigenvector. When A = -1, from 

(1)2 x-2y — 0, giving an eigenvector • 

Considering the eigenvector along the x axis, this implies M is already 
in basic form - since the basic matrices have the axes as their eigen¬ 
vectors. So we test whether M is a direct combination of basic matrices. 

By inspection M = 1) (o - 1) a re ^ ect ^ on m t ^ ie x ax * s followed 

by a shear of factor 2 in the x axis direction (Fig. 6.3). It is worthwhile 

to see how the eigenvector ^ j ^ deals with the situation. Consider a new 

set of axes x lt y x with the eigenvector as the new x x axis. The x x axis 
is thus at 45° to the x axis. The rotation of axes matrix 


u _ / cos 45° sin 45°\ __ 1 / 1 1\ 

sin 45° cos 45°/ J2 \~1 1/ 

is needed to transform the coordinates of the points. 
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Fig. 6.3 

The mapping of the matrix M when referred to the x x y x axes is 
given by RMR -1 (refer to page 57) 

"^(-1 0 (o -1)72(1 1) 

-K-l DC -!) 

-1(11) 

-(1 1)- m - 

Ml is the matrix now referred to the x v y 1 axes. We express this in terms 
of the basic matrices. 

m ‘-(J D(1 - a reflection in the y 1 axis followed by a 

shear of factor - 2 in the x 1 axis direction. The reader is invited to find 
the coordinates of the old unit square with respect to the x v y x axes and 

follow through a reflection in the y x axis and the shear ^ ^ 

verifying that the same transformed figure is obtained. 

(ii) M = (j 2 ). detM= -9. 

The characteristic equation is (2 - A)( - 2 - A) - 5 = 0 from which 
A 2 = 9 and A=±3. 

The equations to give the eigenvectors are 

(2-A)*+y = 0 .... (1) 

5# + (-2-A)y = 0 . . . .(2). 


and 
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When A = 3, from (1) -x+y = 0, giving an eigenvector and, when 
A = -3, from (1) 5#+y = 0, giving an eigenvector • Considering 


the eigenvector the rotation of axes matrix for referring the points 
of the mapping with respect to axes x v y x inclined at 45° to the x> y axes 


■£(-! !)■ 

Thus M referred to the new axes is given by 

!)(! -2)72(1 

-(0 


By inspection 


«.“(o i)(i -?) 

-0. % -?) 


- a reflection in the x x axis, followed by a shear of factor f in the x x 
axis direction and with an enlargement factor of 3 . 

If the eigenvector ^ _ 5 ^ ls considered the rotation of axes matrix is 

1 (\ _5\ 

- 7 = ( - 1 ) since tan a = - 5. The reader is invited to check that the 

n/26 \5 1 / 

same transformation is produced (but with slightly different basic 
matrices) with this eigenvector as the new x x axis direction as was found 

with the eigenvector (as the x 1 axis direction. 


< a ) M -i(" »)• 

We can find the eigenvalues by considering (11 - A)(14 - A)-4 = 0, 
from which A 2 - 25A +150 = 0 and A = 10 or 15. 

The equations to give the eigenvectors are 

(ll-A)*-2y = 0 . . . .( 1 ) 

-2* + (14-A)j/ = 0 .... ( 2 ). 


and 
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When A = 10, from (1), x - 2y = 0, giving an eigenvector 



When A = 15, from (1), - Ax - 2y = 0, giving an eigenvector 



( 


Consdering the eigenvector 



the rotation of axes matrix is 


cos a 
-sin a 


sin a\ 
cos a ) 


where tan a = J; i.e. R = ~j^ 



The transformed matrix of the mapping, when referred to the x 1 y 1 
axes when the direction of x ± with respect to the x axis is tan -1 (J) is 
given by RMR -1 

_ J _/ 2 1 \ 1/11 - 2 \ J _/2 - 1 \ 

“s/n-1 V5\-2 14/ V5 \1 2) 

_ J _/ 20 10\/2 - 1 \ 

“25 V -15 30/ \1 2) 

1 /50 0\ 

— 25 \ 0 75/ 


-(o S) 

- a double stretch of factors 2 and 3 with respect to the x x and y x axes 
respectively. 


6.6 Classification of a mapping 

We now begin to classify the mapping of a matrix by the consideration 
of its eigenvalues and eigenvectors, using exercises for the reader. 


1 

1. Show that the matrix ^ 




gives a mapping in which 


points are reflected in the line 2y — x . (Find the eigenvectors of the 
matrix and show by a suitable rotation of axes that the matrix mapping 
is a reflection in one of the new axes.) 


Verify the determinant of the matrix is - 1. 

Verify also that the following matrices have a determinant of - 1 and 
find what mappings they produce. 
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2. Show that the matrix M = ^ ^ has determinant - 9 (i.e. of the 
form - k 2 ). Verify it has eigenvalues of 9 and - 1 and eigenvectors ^ 
and ^ respectively. By a rotation of axes through 45° correspond¬ 
ing to the eigenvector ^ show the matrix transforms into the matrix 
= when referred to the new axes. Verify this matrix 

M-,^ * s (o l) (o -1) ~ a re ^ ect ^ on t ^ ie new x followed by a 
linear stretch of factor 9 in the new x axis direction. 

A determinant of the form - k 2 indicates the presence of a reflection 
component in the mapping which becomes clear when the matrix has 
been referred to suitable axes. 

Investigate the following matrices which have determinants of the 
form - k 2 . 


(i) 




If we now stop for a moment and look at the results we have obtained, 
we may note that, after suitable rotation of the axes, all the mappings 
considered reduce to combinations of some or all of shear, double stretch 
and reflection; also that the reduced matrices all have a zero in the 
lower left-hand corner. This latter point is worth investigating to see 
what conditions are necessary for it. It happens when we make the x 
axis an eigenvector of the matrix. Considering the general matrix 

^, if the x axis, represented by the vector is an eigenvector 


a b 
c d 


we must have {a -A)l +b . 0 =0 and c . 1 + (d - A) . 0 = 0.The second of 
these gives c = 0. 

Conversely, if c = 0 in any matrix, these equations show that will 
be an eigenvector with eigenvalue A = a. So we have proved that 


c = 0o(l,0) is an eigenvector. 

The matrices so far considered have all had real non-zero eigenvalues. 
We consider next cases in which the eigenvalues are complex, and 
therefore there are no real eigenvectors. In these cases the mapping 
contains a rotation component which must be separated out before the 
rest can be analysed. 
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6.7 Example with complex eigenvalues 

Consider the matrix M = ^ -1) ‘ ^ ias e ig enva ^ ues 1 ± 2 i. 

We therefore try to express M as M X R where R is a rotation matrix 
and M x has real eigenvalues. 

Taking R = ( C . OS a S * n , since M=M,R^M,=MR 1 we have 
G \ - sin a cos a/ 1 1 

jyj _ /3 — 2\ /cos a — sin a 

1 \4 -1/ \sin a cos a 

3 cos a - 2 sin a - 3 sin a - 2 cos a 

4 cos a — sin a — 4 sin a - cos a 

and a has now to be chosen to make the eigenvalues of M x real. 



The condition for this in the general case is that 

/*)-• 

shall have real roots; that is (a -X)(d-X) -bc = 0 has real roots; that is 
A 2 -(a + d) A + (ad -bc)= 0 has real roots, for which the condition is 


(a + d) 2 - 4 (ad - be) ^ 0 . 

This gives a 2 + lad + d 2 - \ad + 4 bc^ 0 , or 

(a-d^ + ibc^O. . . . (4). 


This condition is not difficult to satisfy - in fact, it being an inequality, 
when applied to M x it will give a whole range of values of a. We may 
therefore consider whether a particular choice of a within the possible 
range may be advantageous. For example, if in (4) we have c = 0, the 
inequality is certainly satisfied, and we then have a matrix M x with zero 
in the lower left hand corner which, as we have shown above, (p. 64) 
already has the x axis as an eigenvector. We do not need therefore, to 
perform any change of axes; M x is already in a form in which it can be 
reduced to a combination of shear and linear stretch. The working is as 
follows: 

Choose a so that 4 cos a-sin a = 0; i.e. 

. . 4 

tan a =4, sin a = —j= , cos a 


1 

= Jr7' 
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Then 


1 


8 


Mi ~n/17 C 0 


-12 -2 
-16 


1 

= Vl7 


JV7 

-1 


( o - 17 ) 

1 /5 0\ /I 2-8\ 

7 VO 17/ VO 1 / 


< 


0 


9 


0\/_L Ox/ 1 2-8 

)Ul7 ( 

W 0 Vl7 A 0 


;) 


which is a shear of factor 2*8, a double stretch with factors S/J17 and 
J17, and a 180° rotation, in succession. 


6.8 Matrices with zero eigenvalues 

It remains to consider matrices which have one or more zero eigen¬ 
values. Take, for example, the matrix ^ ^ > which has eigen¬ 

values 0, 5, or ^ j _ 2 ) > which has the eigenvalue 0 only. It is left as 

an investigation for the reader to show that in these cases it is still 
possible to perform a rotation of axes to make the x axis an eigenvector, 
and to remove stretch components as in the examples above. The first 

type then reduces to one of the forms or ^ ^ , according to 

whether the new x axis, belongs to the zero or the non-zero eigenvalue; 
these give the first two projections shown in Fig. 6.4. 

The matrix with eigenvalue of zero only reduces to ^ which 

gives the two-stage projection shown in the third diagram of Fig. 6.4. 
The details are left for the reader to investigate. 



Fig. 6.4 
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6.9 Further investigations 

1. For the reader who is familiar with the ideas of null-space and 
range of a linear mapping: relate these to the eigenvalues and eigen¬ 
vectors of a singular matrix. 




2. As an alternative approach to the analysis of mappings produced 
by 2 x 2 matrices with non-zero determinant consider the unit square 
being mapped into a parallelogram, which can then be reduced to the 
square again by the sequence of operations shown in Fig. 6.5. Investi¬ 
gate whether, given the original matrix, the necessary matrices to per¬ 
form this reduction can be found. As an extension, is this the only or the 
best possible order for the operations? 


3. In the analysis in this section we have required the shear com¬ 
ponent to be in the direction of the x axis. If we allow it to be in any 
direction, can we eliminate the need for a double stretch, making do with 


a homogeneous enlargement 



instead? If so, we can state that 


every non-singular linear mapping of the plane is equivalent, at most, 
to a shear and an enlargement, together with a rotation or a reflection. 
Investigate this. 
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7.1 Introduction 

In this chapter we consider the use of matrix methods in networks 
and communications. The application of matrices in this field is becom¬ 
ing extensive. The range of mathematical technique and difficulty is 
wide, but begins at such a level as to admit of matrix methods being 
taught through this topic in the classroom at the early secondary stage. 
Networks have in fact much to commend them as a teaching approach 
to matrices. The ideas, methods and examples in this chapter have been 
written, in the main, as investigations - for the reader to be actively 
involved. These situations are full of extensions and possibilities and 
once the situation is appreciated then with active participation the reader 
should find much mathematical enjoyment from continuing the situa¬ 
tion in his way with his mathematics. 


7.2 Matrix description of a network 


a b 



Fig. 7.1 


This network can be described by the matrix 


abed 
a j 0 1 1 1\ 

b 1 0 0 0 
c\l 001/ 
d \l 0 1 0/ 


= N 


where the nodes a , b y c , d are written as a row and column and entries 
in the table (matrix) they form are a 0 if the two nodes are not joined 
and a 1 if they are joined. Note the 0 for node a to node a . 

In this representation of a network by a matrix, the matrix is just a 
neat compact store of information, describing the connections within 
the network. 
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The matrix N can be called the ‘node matrix* as it describes which 
nodes are directly connected. 


Exercises 

7.1 Describe, by a suitable matrix the network 



7.2 Draw the network to describe the matrix 


(0 1 I 1\ 
10 11 
110 0 
\l 1 0 0 / 


Immediately, important and fundamental questions can be asked 
about the matrix N - is N symmetrical? Why? What do the sums 
of the individual rows mean? 



4 c 
Fig. 7.3 


In Fig. 7.3 the links between nodes are numbered. The ‘link* matrix; 
the matrix describing which links are directly connected to other links 
is 


12 3 4 
1 /0 1 1 0 \ 
2/1 0 1 l\ 
311 1 0 1/ 

4\0 1 1 0/ 


and the ‘node-link* matrix (which need not be square) is 

12 3 4 
a j\ 1 1 0\ 

b 1 0 0 0 \ 

c\0 1 0 ll 

d\0 0 1 1/ 

(The entry is 1 if the station is connected to the link, 0 otherwise.) 
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Can we construct a matrix to show a one-way street system? - such 
as for the map 


The matrix is 



from 


Fig. 7.4 

abed 
a j 0 1 0 1\ 

b 0 0 0 0 
c\l 1 0 1 / 
rf\l 0 1 0/ 


to 


Note the from y and ‘to’ convention, which here needs to be specified as 
the matrix is now not symmetric. This convention is important and will 
be used in later situations. 

For this matrix, what do the sums of the entries in the rows mean? 
What do the sums of the entries in the columns give? 



c 


Fig. 7.5 


How can a matrix show more than one link between two nodes in a 
network or a link around a node? 


a b 
a/2 1 

6(1 0 

c \1 2 


c 



Exercise 


7.3 


. z 1 0 h 

Draw the network described by the matrix II 0 11. 

VO 2 0/ 
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7.3 Addition of Matrices 

Four towns a , b , c and d are connected as illustrated by the network 



Some new roads are planned to be added to the network. These are 

b 



The new network will be 



In matrix form this situation is 

/0 1 1 0 \ /0 0 0 0 \ 

fo 0 1 0 , ,, 0 0 0 1 

0 1 combined with l Q j Q J 

\0 0 1 0/ \0 0 1 0/ 


gives 


/0 1 1 0\ 
0 0 11 
1101 
\0 0 2 0 / 


This method of combining two matrices is akin to addition and we 
write 

✓0 1 1 0\ /0 0 0 0\ /O 1 1 0\ 

(0 0 1 o\ [0 0 0 l\_/o 0 1 l\ 

11 0 0 1/ \0 1 0 0 ~\l 1 0 ll 

\0 0 1 0/ \0 0 1 0/ \0 0 2 0/ 

where in the addition of two matrices we add corresponding entries. 

The network which when added to any existing network leaves the 
network unchanged will be the matrix of zeros, i.e. no new roads 
added. Such a matrix is the identity matrix for addition. 
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Exercise 

7.4 Add the matrices and draw the corresponding networks for 


(\ 0 lx 

/0 1 

i\ 

l 1 0 

v 0 1 

°) 

Vo 2 0/ 

Vo o 

0/ 


7.4 Multiplication of matrices 

Four towns a , b , c and d are connected by bus as illustrated by the 
network 



Fig. 7.9 

I enjoy walking and have found walks between these towns, with the 
best direction for walking, as 


\ 


->■-- 

ur d 

\ 

Fig. 7.10 


The matrix which describes the bus connection is 


B = from 


abed 
a /0 1 1 0\ 

b 0 0 1 01 
c\l 001 
rf \0 0 1 0 / 


and the matrix for the walks is 


W = from 


abed 
/0 0 0 0 \ 
10 0 1 
0 1 1 
d\ 0 10 0/ 


to 


to 
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I like going by bus from one of the towns and then taking a walk. 
What routes are possible? 

Combining Figs. 7.9 and 7.10 we have the net-work as 



Fig. 7.11 


I can take a bus to b and walk back to a from b so a to a is a possible 
route, a to b is not possible as I can go to c and b by bus but I have no 
walk from c or b to b . The matrix which shows routes between the 
towns when a bus journey is taken followed by a walk is, by considering 
the network of Fig. 7.11 


M = from 


abed 
a i 1 0 1 2\ 

b 0 0 1 l\ 
c\0 1 0 0/ 

d\0 0 1 1/ 


to. 


Can we construct the matrix M from the matrices B and W? This 
situation leads to a method of combining matrices. Can we go from a to 
b ? asks the question - ‘is there a town x , so that one can go from a to x 
by bus and then from x to b by walking? We find which towns we can 
go to by bus from a by looking along the row a of matrix B. We may 
travel to b or c. Now if we are able to walk to b, we look in column b of the 
matrix W and find we can only walk from dtob: i.e. from b or c> there is 
no walk to b. Hence a to b does not have a journey by bus followed by a 
walk, and the 0 appears in the a row, b column of the matrix M. Con¬ 
sider the routes between a and d. The row a in B says we can go from a 
to b 04 c and the column d in W says we can come from b or c to d , i.e. 
there are two possible routes from a to d. We can find if there are routes 
and how many by calculation - using the 0’s and l’s - for there is a 
route from a to d only when a 1 coincides in the row a of B and the 
column d of W. The row a in B is 0 1 1 0 and the column d in 

WisO 
1 
1 

0 . 
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The number of possible routes from a to d can be calculated from 
0.0 +1.1 +1.1 +0.0 =2, taking the corresponding entries in row a of B 
with those of column d in W and adding. This is matrix multiplication. 
Once the situation is appreciated, how much easier it is to obtain M 
by matrix multiplication from B and W than having to sort it out from 
the network! Check BW = M by matrix multiplication. Consider WB 
by matrix multiplication and consider the result by both calculation and 
description (non-commutative law for matrices emerges here). Does 
B 2 = (B . B) have any meaning? (B 2 gives the possible routes between 
towns such that a bus journey is followed by a bus journey.) 

What happens if I take a bus journey and then a walk and then an¬ 
other bus journey? Will matrices and matrix multiplication cope with 
this situation? Try other combinations - such as BWW. What com¬ 
bination of B and W will allow me to get back to my starting town? 
(notice principal diagonal of M). What least combination of B and W 
allows communication between all the towns? 

The power of the matrix application is that the rule of matrix multi¬ 
plication applied in these situations gives the routes possible whatever 
sort of journeys are asked. There is no need to trace separate journeys; 
the matrix multiplication picks up the possible routes and supplies a 
final ‘route data' matrix. Note that a bus journey (B) followed by a walk 
(W) is given by BW (and not WB - the reverse order met in transforma¬ 
tions). This (true) order is due to our convention of writing ‘ from* and 
Ho*: reverse ‘ from\ Ho* and the conventional backward multiplication 
is required. 

With the bus journey matrix of Fig. 7.9 what would be the walks 
such that a bus journey followed by a walk gave the same connections 
as that of the bus journeys? This walk matrix would be the identity 
matrix for multiplication. The bus and walk connection would need to be 



Fig. 7.12 


i.e. walks around the towns themselves. The identity matrix is thus 

/I 0 0 0\ 

[ 0100 ] 

10 0 1 0 / 

\0 0 0 V 

- a matrix with l’s along the principal diagonals and 0’s elsewhere. 
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7.5 When is matrix multiplication possible ? 

To find how the ‘shapes’ of two matrices must be related in order to 
multiply them, consider the following situation. 



I live in Bristol and friends live in London. We have agreed to meet 
in Paris, and I shall travel to Southampton and can take the channel 
steamer to Le Havre or Cherbourg. My friends will go to either Dover 
or Southampton. From Dover the channel steamer goes to Calais or 
Boulogne. 

The matrix (Mj) giving routes from B and L to the Channel ports is 
S D to 

from ^ ^ j 2 ) anc * t ^ ie matr * x (M 2 ) giving routes across the channel 

Le Ch B C 

is D (0 0 1 l) * From t ^ ie French channel ports to Paris the 

matrix M s is 

P 
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The matrix, found from the network giving possible routes from B and 
L to the French ports is 


M 4 = 


Le Ch B C 
B ( 1 1 0 0\ 

L \1 11 \)' 


Again M 4 can be obtained from Mj and M 2 by matrix multiplication and 


S D Le Ch B C Le Ch B C 
B( 1 0\ S (1 1 0 0\_5/T 1 0 0\ 

L \1 1 / D\ 0 0 1 1 / Z-^l 1 1 l)' 


Note the position of S and D and how the matrices are combined. The 
network is simple enough to allow the possible routes from B and L to 
Paris to be traced and the number of routes is given by the matrix 


(M b ) = 



This matrix is obtained by using the matrices M 1? M 2 and M 3 as 



The matrix which describes the routes from Southampton or Dover 
to Bristol or London is 


B L to 


from 


S(1 
D VO 


0 


This is the matrix M x with the rows of written as columns. 


The matrix formed from another matrix by writing the rows as 
columns is known as the transpose of the original matrix. The transpose 
of is denoted by M[. By transposing the matrices M 2 and M 4 , 
show 

M 4 = = M 2 M^ 

(note the reverse order of and M 2 ). In general for any two matrices 
A and B, (AB)' =B'A'. 


7.6 Another situation where matrix multiplication is applicable 

The ‘family tree* of the male side of ‘father* and ‘brother* is as 
shown 
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Fig. 7.14 


Construct the ‘father’ and ‘brother’ matrices. 


F a b c d e f father B a b c d e f 
a /0 0 0 0 0 0\ a /0 0000 0\ 

W1 00000\ b 00 1 00 0\ 

c 1 0 0 0 0 0 1 c 0100001 

d 0 1 0 0 0 0 d 0 0 0 0 1 0 

e \ 0 1 0 0 0 0 e \ 0 0 0 1 0 0 

f \ 001000/ / \0 0 0 0 0 0/ 

Why is the matrix B symmetrical and F not symmetrical? Construct 
the matrix for ‘uncle’. 


nephew 


abode f uncle 
a [0 0 0 0 0 0\ 

b 0 0 0 0 0 0 \ 

c 0 0 0 0 0 0 1 

d 0 0 1 0 0 0 * 

e \ 0 0 1 0 0 0/ 

/ \0 1 0 0 0 0 / 


How can we construct the matrix U given the matrices F and B? 
(U = FB - note how this depends on the ways in which F and U were 
written down. If U' had been put down, then we need a different rule 
than we have devised for multiplying matrices; a row on row operation 
if FB is to equal U'. BF would not even give U', since we have from 
U = FB, then U'= (FB)'=B'F'=BF' (as B is symmetrical). So if F' 
had been written down, the connection is U'=BF' with the accepted 
rule for matrix multiplication. What rule would be needed to connect 
F', Band UP). 

(This example shows that there is not necessarily one rule for matrix 
multiplication. The rule ‘row on column’ is considered the easiest and 
is the usual one to turn up.) 
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7.7 Distances around networks 

A salesman needs to travel through all the towns on the map. What 
route should be taken to travel the least distance ? 


e 



The matrix 

a b c d e f 

a / 0 3 2 4 - 10\ 

b 3 0 5 1 2 - \ 

c 2 5 0 - - 5 l represents the distances between 

d 4 1-02 - I two connected towns. 

4 - 2-20 6 
f \ 10-5-6 0/ 


Make up a table which gives the shortest distance between any pair of 
towns via a third town. 


a 

b 


M = 


c 

d 


f 


a b c d e f 

0 5 8 4 5 7\ 

5 0 5 4 3 8 \ 

8 5 0 6 7 12 | shortest distance between a and b 

4 4 6 0 3 8 is via</( = 5). 

5 3 7 3 0 -/ 

7 8 12 8 - 0/ 


How are the matrices A and M related? To find the shortest distance 
from a to b (via some other town) we take the distance from a to the 
other towns excepting b (top row of matrix A) with the distance to b 
from the other towns (second column of matrix A) and add the corres¬ 
ponding distances - the smallest will be the shortest distance from a to 
b via another town. Top row of A is 0324-10 and second column of A 
is 3 
0 
5 
1 
2 
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The corresponding additions are 3, 3, 7, 5, - so accepting the 

initial two 3’s which represent the distance a to b, the shortest distance 
is 5 (via d ), the other distance 7 is via c and the dashes show it is impos¬ 
sible to get between a and b via other towns. M can then be built up in 
this way. We can define M as being A # A where the operation * implies 
taking the smallest value of the corresponding additions of two elements 
in appropriate row and column entries of A. 

What does A*A*A mean? 

Can the shortest route (if any!) between all the towns be found in this 
way? This situation is similar to critical path analysis. 

Is there any situation in which taking the maximum (rather than 
minimum) of the corresponding additions could be relevant? 


7.8 Incidence matrices 

Figure 7.16 shows four stations a , b , c, d connected by 5 lines labelled 
1, 2, 3,4 and 5. The ‘station-link’ matrix for this network is 



Fig. 7.16 


12 3 4 
a j\ 1 0 0 
b 10 11 
c\0 1 1 1 
d\0 0 0 0 


5 

°\ 

0 ] (Station-link implies stations as rows and 
1 J links as columns.) 

1 / 


This matrix is an incidence matrix for it is made up by considering where 
routes and points are incident. 


Exercises 

7.5 (a) What do the sums of the individual rows of A signify? 

(b) What do the sums of the individual columns mean? 

7.6 How can the closed loops in the network be found from the matrix 
A? 

7.7 (a) What links need to be cut to isolate the station a ? 

(b) What links need to be cut to isolate a and b (together) from 
the network? How can these links be found from the matrix? 
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7.8 From the transpose of A, A'. Compute the matrix AA'. What 
information is given by the entries of AA' ? 


7.9 For the network Fig. 7.17, label the regions a, /?, y; y being the 
whole outside region. Find the ‘link-region’ and ‘station-region* 
matrices, noting them P and Q respectively. Compute and interpret 

PP', P'P, QQ'. 


a 



7.10 Given a network we can always describe it by a matrix, but given 
the matrix - can we draw the corresponding network ? 

Draw the networks corresponding to the matrices 


(i) 12 3 4 

a 11 0 1 1\ 

bl1 1 0 0\ 

c\0 1 1 01 

d\o 0 0 1/ 


(ii) 12 3 4 
a /l 1 0 0\ 

bl0 1 1 0] 

c\l 0 1 0/ 

rf\l 0 0 1/ 


What conditions are necessary on the matrix for a network to be 
drawn? 


7.9 Dominance Matrices 

In the House Championship Football Cup matches of four Houses 
A,B,C and D the results were 

A had beaten B and D 
B had beaten C 
C had beaten A and D 
and D had beaten B. 

These results can be graphed as 


A B 



Fig. 7.18 
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where A->B implies *A beats B y i.e. A is dominant over B. We can also 
record the results in a matrix as 


A B C D 
A j 0 1 0 1\ 

BOO 1 0 

C\1 0 0 1/ 

D\ 0 10 0/ 

where the entry in row z, column j (a {j ) is a 1 if team i beats team/ and a 
0 if team i is beaten by team /. Since a team does not beat itself the 
entries in the leading diagonal will be zeros (i.e. a u = 0). How can we 
find from the matrix how many games a team has won? 

In the results A has beaten B and B has beaten C - so we would 
expect A to beat C. But this is not so - an upset of form has occurred! 
Therefore the relationship is not transitive (if a {j = 1 and a jk — 1 then 
a ik may or may not be 1 - C has beaten A and A has beaten D so C 
should beat D - and this time has!) 

If A beats B , then A is said to be directly superior to B . If A beats 
B and B beats C then A is said to be indirectly superior to C (even if C 
does beat A). 

For the matches the ‘indirectly superior’ matrix M for the houses 
(from the network) is 

A B C D 
A / 0 1 1 0\ 

_ B /1 0 0 11 Note. C is indirectly superior to B twice 

CIO 2 0 1/ (via D and A). 

D\ 0 0 1 0/ 

Is it possible to find the matrix M from the matrix R? In row A of 
R, A has beaten (is directly superior to) B and D , and in column C we 
find C has been beaten by B so A is indirectly superior to C. Also in 
column Bj B has been beaten by A and C so via Z), A is indirectly 
superior to B. Thus by taking the rows of R with the columns of R and 
performing matrix multiplication we can compute the indirectly superior 
matrix, i.e. M ==R 2 . 

The cup is won on points. If a point waS given for each win we can 
find the number of points gained by each team by summing the rows 
of R. From R the points results are 

A 2 

B 1 

C 2 

D 1. 
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Who is to have the cup? If we were to consider ‘indirectly superior’ 
also as a criterion for points, then M( =R 2 ) 

A 2 

B 2 

C 3 

D 1. 

Putting these results together we have a points result as 


A 

B 

C 

D 


4 
3 

5 
2 . 


So C has the cup! 

The sum of the entries in a team’s row of R and R 2 is known as the 
power of the team, i.e. P(A) is the sum of the entries in row A of 
R+R 2 . 

Consider matches between 5 teams resulting in the following directed 
graph. 



B 


Then 


D C 
Fig. 7.19 

A B C D E 


R = 


A /0 
B 
C 


D 


0 111 
1 0 0 0 0 
0 10 0 0 
110 0 


E\ 0 111 



R gives direct superiority or one-stage dominance , R 2 gives indirect 
superiority or two-stage dominance . 

Pts. Powers (R + R 2 ) 


R 2 = 


(0 

3 

2 

1 

°\ 

6 

A 

9 

0 

0 

1 

1 

i\ 

3 

B 

4 

1 

0 

0 

0 

0 

1 

C 

2 

1 

1 

0 

0 

0 

2 

D 

4 


2 

1 

0 

0 / 

4 

E 

7 
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Note we have a top team (A) and a bottom team (C) from the 2-stage 
dominance. 

It is possible to continue this method to find which is the superior of 
B and D ? Has R 3 any meaning? R 3 will give 3-stage dominance. 

/3 3 1 0 0\ 

0 3 2 1 0\ 

R 3 = 0 0 111. 

10111 / 

\2 1 1 1 1 / 

But R 3 says that A is dominant over itself 3 times (via E and B, D and 
B f C and B) and the other teams are also dominant over themselves as 
given by the leading diagonal. As this dominance is meaningless the 
values in the leading diagonal of R 3 will not be counted. R 3 becomes 

Pts. Powers (R + R 2 + R 3 ) 

( 0 3 1 0 0\ 4 13 

0 0 2 1 01 3 7 

0 0 0 1 1 2 4 

10101/3 7 

2 1110/5 12 

We have still not differentiated between teams B and D. If further 

stages are considered more routing through the team itself becomes 
possible. The matrix does not show the teams used, only the number 
of routes between teams and if we do not count for dominance those 
routes through the team itself, we have to find these routes from the 
diagram. Although the matrix does give the number of routes to be 
found, the matrix method is not practical beyond the third stage. 

Exercise 

7.11 Intelligence finds that the communication links between five 
stations is as described in the network, Fig. 7.20. Find the matrix to 
describe the network. If communication is to be destroyed, which 
is the most vital station to destroy? Make sure of your answer by 
drawing the network having eliminated your station. 



Fig. 7.20 
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7.10 Communication Networks 

Communication between four people (each connected at least one-way) 
is described by the graph, 



d 

Fig. 7.21 

The matrix which describes this situation is 

abed 
at 0 1 0 1\ 

a M 0 0 1 l\ 

A== c 1 0 0 1 

d\0 0 1 0/ 

where the entry a {j is a 1 if i can communicate withy and a 0 otherwise, 
i.e. 012 = 1 as a can communicate with b whereas a 21 — 0. Note a u — 0; 
i.e. a person cannot communicate with himself! The matrix A describes 
one-stage communication between two people - and A 2 gives 2-stage 
communication. 

From this example (and others) it can be seen that there is always one 
individual who can communicate with every other individual in 1 or 2 
stages. The situation in which this happens must be such that the 
entries exhibit the properties, 

(i) «« = 0 

(ii) either a ij = t or a H = 1 (there must be one communication link 
between any two persons). To prove that at least one individual 
can communicate with every other individual in 1 or 2 stages it is 
necessary and sufficient to show that in A + A 2 there exists one 
row with all elements different from zero (except the element 

re¬ 
consider the situation in which two individuals i and j cannot com¬ 
municate in 1 or 2 stages. Then 0 ^ = 0 in A and also c ij =0=A i A i 
where A* is the row i in A and A j is the column j in A and thus the 
product of the row and column (for A 2 ) gives zero. Now a u =0 implies 
a H = 1 (from initial conditions). As A t A j = 0 then if the element a ik = 1 
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(the &th element in row i) then we must have the element a kj = 0 (the 
>th element in row k) so that the product on matrix multiplication is 0. 
If a k j = 0 then a jk = 1; i.e. a ik = 1 implies a jk = 1. Hence we have the 
situation that in a matrix A with n persons, the elements are such that 

a 3l + a i2 + a iz + • • • + a jn ^ a i\ + a i2 + a i% + • • • + a in + 1 

since we have shown if a i5 = 0 then a H = 1 and if a ik = 1 then a jk = 1 and 
also one of a i2 is considered zero (hence the 1 on the end) i.e. if one 
individual cannot communicate with another individual in 1 or 2 stages 
then in the matrix A 

a n + a j 2 + a j a + • • • + a jn > a n + a i 2 + n iz 4-... 4- a in (1). 

Now let m be the individual for whom the sum of the elements in the 
mXh. row of A + A 2 has the greatest value of all the rows. Then m can 
communicate with all individuals in 1 or 2 stages, for if he cannot (1) 
above assures us there would exist another row with a greater sum (since 
we began with the premise that one individual could not communicate 
with another individual in 1 or 2 stages and were led to result (1)) and 
this contradicts the statement that the row m has the greatest sum. 

Exercises 

7.12 if a b c d e 

a j 0 1 0 0 1\ 

b 0 0 0 1 1 1 

A= c 1 1 0 0 0 

0 1 0 1/ 

* \0 0 1 0 0/ 

(i) draw the directed network 

(ii) calculate A 2 

(iii) is the network such that a u = 0 and either a i} = 1 or a n — 1 

(iv) which individual can communicate with all others in 1 or 2 
stages, as selected by the ‘greatest sum’ (check result on graph) 

(v) which other individuals (if any) can communicate with all 
others in 1 or 2 stages. 

7.13 Rumours: Mrs. A chatters to Mrs. B, Mrs. C and Mrs. F. Mrs. C 
chatters to Mrs. A , Mrs. E and Mrs. D. Mrs. E chatters to Mrs. C 
and Mrs. F . Mrs. F chatters to Mrs. A and Mrs. E> and Mrs. D 
chatters to Mrs. C and Mrs. E . Draw a ‘chatter diagram’ for these 
communications. Construct a matrix showing rumours at first and 
second hand and decide who would spread a rumour the quickest 
throughout the group. 
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7.11 Sociometric Data 

The social attitude between individuals may be summarised into 
three categories - acceptance, indifference or rejection. From a study 
of the social habits of individuals it is possible to develop a sociogram - 
a diagram which indicates the feelings of one individual towards an¬ 
other. In such a diagram: 

a->b ; a indicates acceptance of b - true friendship 
a — b ; a indicated indifference towards b 
a — \b\ a indicates rejection of b - not on talking terms. 

A sociogram between 3 individuals may be as in 


a 



Fig. 7.22 


a accepts b and c 
b accepts a and rejects c 
c accepts a and indifferent to b . 

By means of matrices it is possible to find a method of quantifying 
such a sociogram to show 

(i) a number status for each individual in the group taking into 
account his choices and rejections and also the kinds of choices 
of others in which he is involved. 

(ii) ranking an individual showing his place in the group giving an 
indication of social adjustment or social rejection 

(iii) determination of cliques within a group 

(iv) establishment of a numerical index for the group so it may be 
compared to other groups, to an optimum value for a ‘perfect' 
group and to compare the same group with itself at a later stage 
when different sociometric information is collected. 

The matrices which describe the sociogram Fig. 7.22 are 

a b c 
a/ 0 1 1\ 

A = acceptance b I 1 0 0 1 

c\ 10 0/ 

where a 1 indicates acceptance and 0 indifference. The rows show what 
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choices a person has made and the columns the persons who have chosen 
the head of the column. 

With A is the rejection matrix 

a b c 

a / 0 0 0\ 

R = rejection b I 0 0 -1 ) 

c\ 0 0 0/ 

where a 0 indicates indifference and -1 rejection. Row and column 
meaning as before. 

Consider 

a b c 
a /2 0 0\ 

A 2 - 6 0 1 1 

f\0 1 1/ 

which gives ‘two-stage’ acceptance, i.e. b is acceptable to c via a . 

A 2 describes relations within a group via another person in the group 
- even showing when a person is mutually acceptable e.g. a (twice) and 
b and c (once). These are the values in the leading diagonal of A 2 and 
are the clues to any cliques. 

Now 

a b c 
a/ 2 1 1\ 

A +A 2 = bl 1 1 l) 
c M 1 V 

gives the ‘one’- and ‘two’-stage relationships between this group of 
individuals. 

The total of the columns give the number of ‘one’- and ‘two’-stage 
inward relationships (i.e. the individual has been chosen by these 
number of people) enjoyed by each individual - a> 4 and b and c> 3. 
The row sums give the outward relationships between the individuals: 
a , 4 and b and c, 3. 

If we now use R and compute 

a b c 
a / 2 1 1\ 

A + A 2 + R = bl 1 1 0) 

c\ 111/ 

we have a matrix from which social values can be determined - the 
column totals giving number of inward relationships and the rows totals 
outward relationships. A social index can be found for each individual 
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by combining his row and column total - and if it is considered that 
outward relationships are more important than inward relationships 
we can multiply the row totals by two - before adding to the column 
total. Here then the social index of an individual is the addition of his 
column total with twice his row total in the matrix A + A 2 + R. In the 
example 

# = 4+ 2(4) = 10 
6 = 3 +2(2)= 7 
c=2+2(3)= 8 

and these are the social indices for each individual. The composite 
social index for the group is 10 + 7 + 8=25. Ideally when each indi¬ 
vidual accepts all others, the group index is 3 n\n -1) when n is the 
number of individuals in the group - i.e. for a threesome 54. 

Cliques can be discovered from mutual choices. If we find the 
individuals who have the greatest social index and from A 2 the indi¬ 
viduals who enjoy most mutual choice, then from A we can decide who 
those individuals are and decide if they form a clique. 

Exercise 

7.14 Investigate the sociogram 



(The methods given in 7.11 were suggested by an article by F. Schippert 
in the American Journal, School Science and Mathematics , Dec. 1966.) 
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SUCCESSIVE EVENTS, STOCHASTIC PROCESSES 
AND MARKOV CHAINS 


8.1 The Use of Matrices in Life cycles 

The following example, although somewhat artificial, shows how 
matrices can describe situations in which events recur. Other more 
realistic situations will be dealt with later, but this is a good situation 
with which to begin. 

A certain species of animal has a life cycle such that f of the animals 
in their 1st year of life survive their first birthday; of these, | survive 
their second birthday and live into their third year; no animal lives 
more than three years. 

During the third year of life, reproduction of the species takes place, 
and during this year three animals (on the average) are born for each 
animal alive. If we start with 3,000 animals, with 1,000 in each year 
group, how will the animal population fluctuate over future years? 

The starting situation working in thousands and writing as a column 


vector is 



for the 1st, 2nd and 3rd years respectively. After the 


first year, the thousand animals in their first year will become f . 1,000 
into the 2nd year (since only f of the animals survive their first 
birthday). The thousand animals in the 2nd year become J. 1,000 
in the third year and the 1,000 animals in the third year reproduce 
3 x 1,000 to go into the 1st year and no animal lives into the 4th year and 
beyond. 


The situation now is 



for the 1st, 2nd and 3rd years. 


After the second year in the life cycle, the situation is 



These values are found from the first year values by a calculation 
explained below, which can be shown schematically in the diagram. 

89 


D 
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We have taken the first year value, multiplied it by f and placed it in 
the second year, the 2nd year value is multiplied by \ and put in the 
3rd year, and the 3rd year value is multiplied by 3 and put into the 1st 
year. A matrix has the property of multiplying values and transposing 
the places when multiplied with another matrix and it is therefore 
possible to represent the life cycle of this species of animal by a matrix. 


( 0 0 c\ /x\ f cz \ 

a ^ ^ )( )=(##) - note transposition of x , 

0 b 0 \byJ 

and z and the positions of a, b and c in the matrix. If the matrix ( y 

\s- 

® /0 0 c 
and the matrix ( a 0 0 

VO b 0- 

/0 0 3\ 

as ( f 0 0 ) with the proportions surviving in the year groups, then 


we can find the number of animals in the year groups of the second 
year by 


{0 0 3 

t o o)(t: 


I o 

/0 0 3 


The matrix M 


the animals. 


/U U 

■ I 0 0) 

Vo i o/ 


completely describes the life cycle of 


After three years the situation is given by MI 2 ) that is 



We have arrived at the starting situation and the cycle then repeats. 
Show M 3 = I. 
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Examples 

8.1 Find the total number of animals in each of the three years of the 
above cycle. 

8.2 Investigate and describe the life cycle represented by the matrices. 



is 


8.2 Survival or Extinction of the Species 

A life cycle is represented by the matrix ^ ; that is the species ii 

such that i survive their 1st birthday and live into their second year and 
for each individual alive in the second year 1 is born and none of the 
species lives longer than 2 years. How does the population fluctuate 
over future years? 

Assume that 100 of the species is alive in the 1st and 2nd years as a 
starting situation. 

Total 


0 1 
0 


“li 

4 


=T 1 

4 , 


After 1 year the position is ^ 

After 2 years the position is 

After 3 years the position is ^ (^j ~ (^-) 

After 4 years the position is ^ ^ = {^j 

Will the cycle repeat? 

After 2 years the position can be seen as 

6 J)(K o)’(0- 

After 3 years the position is ^ ^ (l) 

After n years the position is ^ ^ ^ . 


5 

64 
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From the pattern of the above results the total population of the 
species as given by ^ can be seen to be decreasing as n 


increases. 


Example 

8.3 Investigate the situation as described by the matrix ^2 • Can 

you conjecture as to the future of a species in the general case as 
described by the matrix ^ ^ for various values of a and b ? 


8.3 ‘Steady State’ Populations 

We investigate the yearly population of a species as described by the 


matrix 


“(I o)- 


The species is such that for every three individuals in 


their first year one individual is born, (represented by the | in the 
matrix) and ^ of those in their first year survive their 1st birthday and 
live into the second year. During the second year of life each individual 
bears two offspring, and no individual survives the second year. 

Assuming 100 individuals in each of the 1st and 2nd years as a starting 
situation, after 1 year the position is 



since for every 3 individuals in the first year 1 is born to go into the 
next first year and | survive into the second year of life, and of those in 
their second year, 2 individuals are born for each in their second year 
and none survive into the third year. 

The matrix gives the first year situation as 



corresponding to 2f (hundreds). 

For the second year the situation is 
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Or from the matrix 


(i imtrMt) 


For the third year the position is 


(l o)(|)-(||) tota12 ^ 


And for the fourth year the position is 


(I ?)(«)=(¥) <-*■ 


Collected together the number of each of the species in each year of 
its two year cycle and the population in each successive year is 


1st year of life 

1 


1 = 2*33 


2nd year of life 

1 

i=0*33 

1=0*78 


total in year (100’s) 
2 

2f = 2*67 
2f=2*22 




fx=0-68 


2ff = 2*41 


There is much fluctuation in the values here. It appears that the 
population is neither decreasing or increasing and may be tending 
towards a limiting value - ‘a steady state* value. Can we find this value? 

If the population of the species does reach a ‘steady-state* after some 
period of time, then the situation would be that if the numbers in each 

year at this stage were then ^ (^j = (^j that is, the ‘transition* 
matrix ^ describing the change in population from year to year 
has not altered the value of the vector - the ‘steady state*. 

This vector must be an eigenvector of the matrix. In this situation we 
see the importance of the eigenvector; it gives those values of the situa¬ 
tion, wherein after the matrix is continually applied, the matrix does 
not alter the eigenvector values. 


If T is the matrix and x the vector then the steady-state 

position is Tx = x where x is the eigenvector of the matrix T. 


k (! o) 
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If A is an eigenvalue of the matrix T = 


f the matrix T = ^ ^ 1 

(l o)GH;)- 


from which 


^x + 2y = Ax 
i# + 0 .y=Ay 


(^-A)# + 2y = 0 
^# + (0 - A)y = 0 


For these equations to have non-zero solutions we must have 


that is -A(^-A) -f = 0; A 2 -^A-f = 0 and A = 1 or = -f. 

From (1), when A = 1 the eigenvector is given by - l)x + 2y = 0, i.e. 
§# = 2y or x = 3 y and as a vector by . 


Bu t since (J q)( 3 *)=( 3 *) 


we can only deduce that when the steady 


state is reached the number of the species in the first year of life is three 
times the number in the second year. We are unable to find the actual 
values of the population in their first and second year of life when the 
steady state is reached. 

To be able to find these values of the steady state we need to be able 
to compute successive powers of the matrix, which will give the popula¬ 
tion values for each successive year, and then endeavour to deduce from 


these values the steady state numbers. Since T = ^ and the popu¬ 

lation for the first year is given by o) (l) an< ^ t ^ ie secon( * y ear by 
^ and the nth year by ^ ^ we require a method of 
computing f1 Q J . 
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8.4 Power of a Matrix 

The other eigenvalue of T was -f and the corresponding eigen¬ 
vector is given by (i + §) x + 2y = 0, that is x= -2y> or by the vector 


check fi \ 




The two eigenvectors of T are given by the vectors and ^ 
Construct the matrix P by writing the eigenvectors as columns , that is 

'-(? D- 


Then the inverse of P is P -1 = ^ 


and computing P _1 TP gives 



This is a matrix of diagonal form with the eigenvalues of T as the 
entries on the major diagonal. (It will be proved later that using the 
eigenvectors of a matrix in the above way will always give a result of 
this form.) 


Let D = (q _ |) then P _1 TP = D. 

We require to find T n . 

Now P(P~ 1 TP)P~ 1 = PDP -1 , pre- and post-multiplying the equation 
P _1 TP = D by P and P -1 respectively, that is 

(PP- 1 )T(PP- 1 ) =PDP 1 
(by the associative law) and thus 

T = PDP -1 . 

Then 

T 2 = (PDP 1 )(PDP 1 ) = PD(P 1 P)DP 1 = PD 2 P _1 . 

T 3 = PD 2 P -1 (PDP -1 ) - PD 3 P-! 
and hence (by induction) 


T n = pDnp-i 
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It is easy to compute powers of a matrix in diagonal form since 


Co -iMi -86 

(o -f) = (o (-|) 3 ) 


4)~(o (-°i) 2 ) 


and by induction 


(J 4)■-(■ 


l o 

o (-!)" 


Hence 


T "-(i i)(o (-fr)K-i 3 ) 

_/3 -2(-fJ-W 1 2\_ 1 /3+2(-f r 6-6(-|)"N 

v (-f)v 5 \-1 3;~ 5 \i-(- f)» 2+3(-t r)- 

Now as « (the number of years in the life of the species) increases from 
the starting situation ^ , T n tends to the limiting state 

i(j 2 ) ( 1 ) (since( -§) n -*0 as n->co) 


Thus the population will tend towards a limiting value of § +f = 2*4 
(hundreds), distributed in the proportions (f, f) = (l*8, 0*6) in the first 
and second years of life. 

Compare this result with the calculated values for the first four years 
of the life cycle, and note the values (f, f) agree with the eigenvector 
of the matrix. In this example the eigenvalue A = 1 and a steady state 
exists; it is not difficult to deduce that if A>1 the population will be 
increasing and that if A< 1 the population is decreasing. 


Examples 

8.4 If X ± and A 2 are the eigenvalues of the characteristic equation of the 
matrix T = ^ ^ with and the corresponding eigen¬ 
vectors, then the matrix P constructed by writing the eigenvectors 

as columns is ( Cl . Show P~ 1 TP = 

\^2 ^ 2 / 

This establishes, as a general result, the method used in the 
example to reduce T to diagonal form. 
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8.5 Let D = P^TP then T n =PD^P 1 



Insert the matrix P and compute the full result of T M . Hence 
deduce that if A x = 1 and A 2 < 1 a ‘steady-state* is eventually reached; 
if Aj> 1 the population will continually increase and if A X <1 and 
A 2 < 1 the population will continually decrease and become extinct. 


8.6 In the following examples describe the life cycle as given by the 
matrix and then consider the future population expectancy of the 
species. If a ‘steady-state* value is to be reached, find this value. 


(a) 




(Investigate how these matrices are devised - the general matrix 
is of the form ). 


8.7 The theory of 3 x 3 matrices describing life-cycles is more com¬ 
plicated. The reader is invited to study matrices of the form 

(i) /0 0 c\ and (ii) /0 d c\ 

[a 00 (a 00; 

Vo b 0/ Vo b 0/ 

e.g. 

/0 0 4\ /0 1 6\ 

0 0) and (j 0 0 J. 

Vo i 0/ Vo i 0/ 

These will be found to describe life cycles in which there will be 
only one stable distribution of proportions of the population in each 
year of its life and there is no total population which the species 
will finally tend towards. (Since the only eigenvector comes from 
A = 1 and it is impossible then to diagonalise the matrix and 
consider powers as «-»oo). 


8.5 Stochastic Processes and Markov Chains 

In the previous section we have seen how a matrix may describe the 
outcome of successive events (the reproduction and dying in the life 
cycle of a species). The matrix is often called the transition matrix as it 
describes the situation after each operation, A sequence of events where 
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the outcome is dependent on some chance element is known as a 
stochastic process . In these events probability is involved and the entries 
in the transition matrix will be probabilities. We examine some examples 
of stochastic processes. 


(i) A black bag (#) contains 2 white and 2 red balls and another black 
bag (y) contains 3 white and 1 red ball. If both bags are held out, what 
is the probability a white ball is drawn out of one bag? There are two 
events to be considered here; one, the choosing of bag x or y and second 
the choosing of a ball in the bag. The probability matrix for choosing a 

bag is^ since it is assumed that each bag has an equal likelihood of 

x y 

being chosen. The probability matrix for the balls is^ If \ \ since in 

x the probabilities of choosing a white or a red is \ whilst in a bag y the 
probability of choosing a white is \ and the red, J. By combining these 
matrices we shall have the probabilities of choosing a white or red ball as 


U i)\V Uxi+ixi/ \t/ 


Hence the probability of choosing a white ball is 5 in 8 times. (Note it is 
assumed also in this situation that if the experiment is continually 
performed, to see how it compares in practice with the theory, the ball 
chosen is replaced before the next draw is made.) The matrix multi¬ 
plication gives the probable outcomes since as the probability of choos¬ 
ing bag x is J and in bag x there is a probability of f of choosing a white 
ball, so there is a probability of J in this case of getting a white ball and 
for bagy the probability of choosing it is \ together with the probability 
of £ of obtaining a white ball having chosen y, that is probability of f 
of getting a white ball from y; the total probability of obtaining a white 
ball is i+f =f. A ‘tree diagram’ shows how the probabilities can be 
obtained and how matrix multiplication performs the same operation. 


AW—W 

x i - *iR 

AW —-f W 
y\^\R -*£R 

The matrix multiplication deals with the situation since we have two 
successive events of which the transition probabilities are known. 


(ii) From a vending machine I can get nut, raisin or plain chocolate 
bars. If I had a nut bar I am equally likely to choose a raisin or plain the 
next time. If I had the raisin I am twice as likely to choose the plain 
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as the nut for the next; and if I had the plain, next time I would always 
get the nut bar! What is the matrix which describes my choices ? 

N R P 1st choice 
2nd N /0 i 1\ 
choice R ( \ 0 0 ) 

PH to/ 

If my first choice is a nut bar, what is the probability I have a nut bar on 
my 3rd choice and what is the probability I have a raisin bar on my 4th 
choice? Note that the starting situation is given and also how the matrix 
multiplication deals with the probability outcomes and check it is 
correct by general probability methods. 


1st choice: probability vector 


W* 

W p 


2nd choice: probability vector = 


3rd choice: probability vector = 



3rd 4th 



Nut = i + i = |; Raisin = Y2 + \ = i- 


(iii) Today is the first day of term. A maths and physics student 
knows that today there is an even chance of working at maths or physics. 
But tomorrow will be different! If he works at maths today there is only 
a probability of ■§■ that he will work at maths the next day (and therefore 
a probability of f that he will switch and do physics). On the other hand, 
if today he works at physics there is a probability of f that he will do 
physics the next day. These probabilities continue for subsequent days. 

At the beginning of th^ 1st day the probability matrix is 

1st day 
M 
P 
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At the beginning of the 2nd day the transition matrix is 

M P today 

M (- -\ 

tomorrow p ( £ %] 

(where the entry f states that if the student did physics today there is a 
probability of f he will do maths tomorrow). 

What is the probability of doing maths on the second day? 


(I !)(IMS)' 


Again check the matrix multiplication deals correctly with the pro¬ 
babilities 


1st 


M(iy s 


P(i)' 


2nd 

(ro) 

^(!) 

-M(\) 

-pm 


What are the probabilities for the 3rd day? 

(I {)’(!)=(} {)©=(!)■ 

Questions such as (ii) and (iii) are answered by taking powers of the 
matrix. In example (iii) what is the probability of the student working 
at maths at the end of the term - after say 100 days? To multiply the 

itself 100 times we use our earlier method of comput- 


matrix 


|) b yi 


ing powers of a matrix. 


If A- 


■(}}) 


or 


i.e. 


and 


so 


then the eigenvalues of A are given by 
(i-A)(f-A)-A=o 

5A 2 -4A-1=0 
(5A + 1)(A-1)=0, 

A = 1 or 
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For the eigenvectors, from (f-A)*+fy=0 when A = 1 the eigenvector 
is 2x -y = 0, or the vector • When A = - £ the eigenvector isx+y = 0, 

or the vector ^ . 

Let P be the matrix whose columns are the eigenvectors, that is 
P = | . Then from work in the previous section 

p - ,Ap - D -(J 4) 

which is a diagonal matrix with the eigenvalues as principal diagonal 
entries. As before A n = PD "? -1 


-(2 -1)(0 (HO-) 1 (2 
_ 3 \2 -Mm2 -V 

3 \2 — 2 ( — 2 +(-*)»; 

i-iM)”\ 
~Vt-fM)“ l+K-iW’ 



the 


When n —100 the entries closely approximate to (and if 

term is endless they do!). 

On the last day of term the probabilities of the student working on 
maths or physics are given by ^ = (^tj s0 there is a 2 to 1 


chance he will be doing physics. 

Such problems as these are known as ‘Markov chains* - where the 
probability at a given instance is a function of the outcome immediately 
preceding. 


Examples 

8,8 It has been found that each year f of the population in a town moves 
into the country whilst \ of the country dwellers move into the 
town. Consider the future balance of population between town 
and country, assuming these proportions of change continue. 
Consider 7,000 town dwellers and 7,000 country dwellers as a 
starting position. 
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8.9 On the 1st day of term a lecturer has a probability of of being 
late for his lecture. But if he is late there is a probability of f that 
he will be early for the next lecture, and if he arrives early there 
is a probability of f he will be early next time. If these probabilities 
continue for subsequent days, what is the probability of the 
lecturer being late at the end of term lecture? 

8.10 Consider the general theory behind these examples. Each 
transition matrix is such that the sum of the elements in each 
column is 1 - why? When a transition matrix is compounded with 
another the sum of each column remains 1 - why? 

The 2x2 matrices in the examples must be of the form 

where 0<tf<l and 0<6<1. Show the eigenvalues of this matrix are 1 
and a + b- 1. Since an eigenvalue is always 1, and a+b- 1< 1, each 
situation will approach some stable state (it will not be possible to have 
a situation where one aspect is absorbed by the other). 


(a 1 

Vi-« 


8.6 The use of the eigenvector in stochastic processes 

As with the life cycles we are looking for the stable state when the 
powers of the matrix no longer change the matrix, that is we have reached 
the state when TP=P where T is the transition matrix and P = T W . 
Again we are considering the eigenvectors of the matrix; for the matrix 
P must have entries of the eigenvectors of T so that TP = P. Finding the 
eigenvectors of T will supply the steady-state proportions and there 
will be no need to compute powers of the matrix. 

Consider the worked example in the text of the student studying 


maths and physics. The transition matrix T 


- 81 ) 


vector ofTbe^j. Then we require the situation 


({t)CH3 whe " 


Let an eigen- 


so that 


ix+iy=x and + | y=y 


-sX+%y = 0 and |*-|j = 0. 


Both these equations reduce to the one equation 2x -y = 0. But we 
know that the eigenvector is a probability vector so that x+y = 1. 
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Hence from 2x -y = 0 and x +y = 1, x , y =f. For the other eigen¬ 

vector we know A = - J (check a+b -1 =■§■ +1- -1 = - -g) then 


Hence 

and 


(i !)(;)=-*(;)' 


5 *+ 5 y= - 6 * 


I x + iy = 


-iy. 


and 


|* + fj = 0 

5 x +iy=o. 

Both these equations reduce to x +y = 0. But since (*^ is a probability 
vector, x +y — 1 - so there is no solution for x and y between these 
equations and the only suitable vector is ^ . The stable state is given 
when A = 1 and then 

(i t\/i iW* *\ 

U *A* fMf i) 

(check with earlier result). 

Example 

8.11 Show from the general 2 x 2 matrix ^ ^ ^ ^ that the stable- 
state eigenvector (A = 1) is 

i^±\ 

ci -\-b — 2 
a - 1 

\a + b-Zj 

Consider the examples in the text by this easier method of finding the 
stable-state matrix which the situation will tend towards. 


8.7 Three outcome situations 

Situations giving 3x3 transition matrices may also be investigated in 
this way. In example (ii) (page 98) concerning choice of the nut, raisin 
and plain bars of chocolate, what is the probability of having a nut bar 
on the last day of term? 
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The transition matrix is 

N R P 1st choice 

N/ 0 * lx 

2nd choice R ( \ 0 0 ) 

PH to/ 

and are seeking the ‘steady-state* eigenvector (again assuming the 


a 


term is endless). If I y J is this eigenvector then I \ 

V*/ H 

and dhox+y+z = l. 


l. 

3 

0 0 


X' 

yj=\y 

z/ 'Z- 


We have the equations 

JL 


■%y+z=x 

h=y 

%x + %y=z 
x+y+z^\ 


• (i) 

• ( 2 ) 

• ( 3 ) 

• (+)• 


From (2) # =2y and then in (1) 3z = 5y (which is the same equation as is 
obtained by putting x = 2y in (3)). 

Then from (4) 
and 


2 y +y + 5y/3 =* 1 




5 

14 * 


Hence x=z Ti an( i 

At the end of term the ‘steady-state’ matrix is 


, 6 

6 

6 > 

14 

14 

14 

3 

3 

3 

14 

14 

14 

5 

5 

5 

"14 

14 

14 ' 


and as a nut bar was chosen at the beginning of term, the end of term 
probabilities are given by 

oWi)* 

W \&Jp 

that is a probability of f of ending the term with a nut bar. These also 
show the proportion of nut, raisin and plain bars of chocolate I have 
had during the term. 


, 6 

6 

6 , 

14 

14 

14 

3 

3 

3 

14 

14 

14 

5 

5 

5 

"14 

14 

14 ' 


Examples 

8.12 A lecturer is either late, early or on-time for his lectures. If he is 
late for a lecture, he is always early for the next. However, if he is 
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early it is twice as likely he will be late for the next lecture as he 
will be on-time. If he is on-time then it is twice as likely he will be 
late as early for the next lecture. If he is early for his first lecture 
of term what are the probabilities of his arrival states for the end 
of term lecture? 

8.13 The diagram shows four compartments with doors leading from 
one to another. A mouse in any compartment is equally likely to 
pass through each of the doors of the compartment in going to 
the next. Find the transition matrix and the probabilities of the 
mouse being in the various compartments in the long run. 


— 

2 

-— 

1 

. 

~3~j 

4 


Fig. 8.1 


8.8 Absorbing state situations 

The examples looked at have been of chain situations in which the 
successive application of the transition matrix has shown the outcome 
of the chain at any stage. Such situations are examples of ergodic chains 
or Markov chains . Since the chains continue for as long as we please 
they are chains with non-absorbing states . An absorbing state is a situa¬ 
tion in which nothing further can happen and the chain ends if that 
state is obtained. Examples with absorbing states are more difficult to 
deal with but are often more interesting. A situation of a Markov 
chain with absorbing states is the following. 

A ‘drunk* lives on a house-boat and has to negotiate the gang-plank 
each evening. One pace to the left or right gets him in the water (the 
absorbing state!!). He wants to get home but the probability (after any 
pace) that he takes a pace forward is while there are probabilities of ^ 
that he will take a pace backwards, J he will stagger left into the water 
and \ he staggers right. (An example of this kind is known as a Random 
Walk!) 

The transition matrix for one step to the next is given by 

F B L R from 

FA i 0 0\ 

P= to B\\ \ 0 0 

Mi i 1 0/ 

R\i i 0 1/ 
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The positions L and R are absorbing states, since once in them it is 
impossible to leave them whilst having paced forwards or backwards it 
is possible to move into any of the other states. 

The drunk is just on the gang-plank so the situation is 



for the first step. For the second step: 


m —i 

— HI) 

R(l) 1 

— ~ R (l) 

i 



i —*(A) 

6 

i—m 

1 

/4' 

—m*) 

m)~~ 

i —*(A) 

6 

~BU e ). 


So the probability of 2 steps forward is ^ +ygr=^. Note that once in L 
or R that state continues. The probabilities for the second step come 

A . 

from using the starting probability vector I jf 1 with the transition 


matrix as 



Again the matrix easily picks up the probabilities. Note that if the drunk 
went into the water on the first step then the chain finishes - but we 
assume it to continue - letting the drunk bob up and down in the water 
for each event considered! Again we need to take powers of the transi¬ 
tion matrix to be able to describe the chain as it continues. The form of 
P allows this arithmetic to be relatively easy. Since P has absorbing 
states (the 1 and 0 columns) it can always be partitioned into canonical 
form as 
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r s 


r 

s 


( 


Q__o\ 

R I / 


Where I is an s x s identity matrix, O is an r x r zero matrix, R and 
s x r matrix and Q an r x r matrix. Is it possible to treat the matrix as 

and multiply it by itself just as we do ‘true* matrices? Does 

/Q OWQ OX /Q.Q + O.R Q.O + O.I\ / Q 2 0\ ? 

\R 1/ \R 1) V R -Q+I - R R.O+I 2 / V R + R Q P/‘ 


Q O 
R I 


The reader is invited to verify that the partitioning of matrices and 
multiplying them matrix-wise as blocks does work. (Use two 3x3 
matrices as 



and obtain the result.) 
Then partition them as 


and let 


c > 


b 

! cS 

\fj k 

1 l \ 

e 

\j_ 

1 1 m n 

\A 

h 

' 

l\p q 

' rj 

B; 

( 8 

h)= C; 

(0-: 


0 


=F; (p q)= G and (r)=H. 


***** (c d)(g h) as matrices and verify, when the capitals 

are reconstituted as matrices and multiplied out (care with the right 
order of the capitals) that the results are the same. 


/Q oy _/ Q 2 o\ 

\R IJ “V R + R Q V’ 

/Q oy_/Q 2 o\/Q o\ 

Vr v ~\ R + R Q i/ V R i/ 

_/ Q* OX 

“V R + R Q + R Q 2 1} 
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and by induction 

/q oy / Q n i\ 

l) ~ \R+RQ+RQ* +... +RQ”- 1 Oj' 

Also 

R+RQ+RQ 2 +... + RQ n_1 =R(I + Q + Q 2 +... + Q w_1 ) 


In our example R = 




If we allow the drunk to be able to take an infinity of steps Q w -> 0 as 
w->oo. (Check this by finding Q n - the eigenvalues of Q are 0 and \ and 


hence Q n = P _1 D W P and D n will be 



which 


0 as w->oo.) 


The matrix Q describes the non-absorbing states so that entries in 
Q n describe the probability of being in any one of these non-absorbing 
states after n steps. Hence I + Q + Q 2 + Q 3 + ... gives the sum of the 
probabilities of being in a non-absorbing state as n increases. 


As 0 as oo, 

(Q 0\n / o 0\ 

\R V IR(I-Q )- 1 V 


as n - 


-CO 


since 


as n 


- Q -(J ?)-(! !M-f 1) “ d t) 


Now 

I 

and hence 


Rd-Q)- 1 ^ |)2(* })-(* J). 

After an infinity of paces the transition matrix becomes 


F B L R 

F (0 0 0 0\ 

BO 0 0 01 

i 1 0 

R\t i 0 1 / 

showing that the probability of absorption in each of the absorbing 
states L and R is each J starting from either of the non-absorbing states. 
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8.9 A Game example 

Many games involve probability and absorbing states - winning or 
losing - and such situations can be considered by matrix methods. 

A simplified game of snakes and ladders could use such a board as 



Fig. 8.2 


with rules of up ladder and down snake; no throw of a 6 on the dice to 
start; if too large a number is thrown to land on 6, ‘home*, the player 
remains on that square and once ‘home* on 6 the player remains there 
(the absorbing state). 

We can find the transition matrix of probabilities of moving from 
one square to another square, as 


M = to 


1 

2 

3 

4 

5 

6 

1 

6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

3 

0 

2 

3 

i 

6 

0 

0 

1 

3 

0 

1 

6 

2 

3 

0 

0 

0 

0 

0 

0 

0 

0 

1 

6 

0 

1 

6 

1 

6 

0 

1 


from 


where entries have been worked out, for example: being on square 3 
and remaining on 3 (column 3, row 3), that is a throw with the dice keeps 
the counter on 3; as a throw of 1, counter goes to 4; a throw of 2, counter 
on 3 (down snake); throw of 3, counter on 6; throws of 4, 5, 6 are too 
high and counter remains on 3. So 4 throws from 6 keeps the counter 
on 3 and then the probability of being on 3 and remaining on 3 is f=f. 
Also we see the probability of being on 3 and going to 4 is to 5, 0 
and to 6, One does not go back to squares 2 or 1. 

The probabilities of being on any one of the 6 squares after the first 
throw are given by 



= Mx = 



square 
1 
2 

3 

4 

5 

6 


so the probability of being on square 6 and home after one throw is 
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The probabilities after the second throw of the dice are given by 


M 2 = to 



1 

2 

3 

4 

5 

6 

1 

/(i) 2 

0 

0 

0 

0 

0 

2 

/ o 

0 

0 

0 

0 

0 

3 

f i 

0 

17 

36 

2 

9 

0 

0 

4 

1 * 

0 

2 

9 

17 

36 

0 

0 

5 


0 

0 

0 

0 

0 

6 

\tt 

0 

11 

36 

11 

36 

0 

1 


u u u u / 

11 11 0 1 / 

36 36 u 1 / 


from 


The probabilities in M 2 can be read off, for example, as the probability 
of being on 3 (column 3) and remaining on 3 (row 3) after two throws of 
the dice is ; the probability of being on 3 and moving to 4 after two 
throws is f. It is a useful exercise at this stage to verify these probabilities 
by the ‘tree* method used earlier - then the usefulness of the matrix 
method becomes apparent; it is a direct and explicit routine method. 

The probability of being on any square after two throws is given by 


M 2 x = 



and the probability of being home after two throws is ^. 


This probability includes that of being home after 1 throw plus pro¬ 
bability of being home after the next throw. To find the probability of 
being home in precisely two throws we subtract the probability of being 
home in 1 throw (•§) from The probability of being home in 2 throws 
exactly is given by the entry at the bottom of the column vector 
(M 2 -M)x. Continuing; the probability of being home in exactly 3 
throws is the entry at the bottom of the column vector (M 3 - M 2 )x etc. 

The average length of a game is considered as 
1 • Pi + 2 . p 2 + 3 . /> 3 +... 

where p n is the probability of being home in exactly n throws. 


For this situation the average length of the game is given by 


1 .Mx + 2. (M 2 -M)x + 3(M 3 -M 2 )x + ...+»(M w -M w - 1 )x . (1) 

when w->oo, considering the entry at the bottom of the column vector. 
Now M can be partitioned as 



1 

3 

4 

6 

1 

(i 

0 

o 

°\ 

3 

i 

2 

3 

i 

°l 

4 

i 

I- 3 - 

1 

6 

2 

3 1 

° 

6 

l* 

1 

6 

1 1 
6 j 

’/ 
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0 O' 


/ 6 \ /Q 0 \ 

and we have Q as I ^ § £ j and R as ^ ■§), so M = l^ jJ . As 


3 6 3 

before as we take powers of M, from page (108), 

Q* 0\ 


Thus 


M H n (i-Q n ) T ] ■ 
' d-Q) 11 


M" - M”' 1 = 


Q" 0\ 

, (i-Q") , 
(i-Q) 

Q" _ Qn -1 

-R 


Qn -1 

' (I-Q" -1 ) 
(i-Q) 
o\ 


°) 

ii 


(Qn _ Q«-l) Q 


\i-Q 

i o 

=(Q”-Q"- l )l -R 


\I-Q 


Oj 


Substituting for values of n from 1 to n in expression (1) leaving the x 
as a factor, the average length of a game is given by 

(l. (Q-I)+2(Q2-Q) + 3(Q^ — R \x 

i o Vq ° 

= ('_I_Q_Q2_Q3..._Qn- 1 + w Qn-l\ R L 

V \i^Q / 

o) 1 - 


It can be shown for the matrix Q that Q n -> 0 as »->oo and also 
«Q n -> 0, so the average length of game becomes 


/ I o\ 

r l 

°\ 

( -I\/ \ 

fi-Q 


(i - q) U °r 

1 R 


\(i - Q) a 

°l 
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/i o <x 

Q =(i f i> 

x 3 6 3' 

/f 0 <x 

(¥ 4 2). 

V# 2 4/ 

R /f 0 Owf 0 Ox 

*»(? ♦ >)(x ♦ >) 

/If 0 Ox 

i)(w 20 20 ) = (6 ^ 

\W 20 20/ 


so i-Q=[ i 

JL JL 

3 “6 


and (I - Q) -1 becomes 


_/l i 
~\6 6 


4 0 \ 
6 )• 


We thus have -—— as a 3 x 3 matrix and the average length of a game 


as 



i 0 .) 

i -Q 

, 40 40 


6 — — 

6 6 

; / 


and the entry required is 6 - the average length of the game. 

With the previous example of the drunkard’s walk, we can find the 
average number of paces before he falls into the absorbing states. Here 


R (F ! Q)') sives 0 0 


so the drunkard, on the average, takes only 


two paces before getting wet. 

As an exercise on the work covered here; find the average length of 
path of the snakes and ladders game when the snake is removed; when 
the ladder is removed; when the snake and ladder are removed. The 
answer in each case is 6, since the 6-sided dice can always finish the 
game from any position on the 6-square board. But consider the situa¬ 
tion on a 9-square board. 

The situations given here are intended to indicate the scope of matrix 
methods. There are other means of dealing with these situations and 
obtaining the answers found, but the application of matrices describes 
the situation neatly and readily produces the answers. The technique 
of matrix multiplication turns up surprisingly often. 





9 
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9.1 Quadratic forms in terms of matrices 

If x = (^j then x' = (x y ) and x'x = (x y) (^j = x 2 +y 2 , that is the 
circle x 2 +y 2 = a 2 can be represented in matrix form by 

(* j 1 ) (*)-«* 

»(J 50-* 

Similarly the equation x' ^ 3 ) x = ^ re( ^ uces as follows 

<* ^(3 5)0- 1+ • ■ • •<')■ 

<* >'l (ivti) - 14 

6 a ; 2 + xy + 3 xy 4- 3jy 2 = 14 

6 a ; 2 + 4a;3;4 3 ^ 2 = 14 . . . . ( 2 ). 


and 


The matrix notation used in (1) for the quadratic equation ( 2 ) shows 
another use of matrices. Scalars, vectors and matrices can be regarded 
as successively more complicated types of number, each having a 
variety of applications. Just as ordinary (scalar) numbers can be used to 
describe the size of a class or the temperature of a room, and a vector to 
describe the position of a point, or the direction of a line, so a matrix 
may be used to describe a rotation of axes or a quadratic form. The use 
of the matrix to describe a quadratic form has wide application in 
geometry and its importance lies in its easy extension into 3 or n 
dimensions. 


Examples 

9.1 What is the equation represented by x'Ax = 1 when 


(3 3)= (e 3 ) 5 (2 3)* 

113 
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9.2 Can you find another matrix A to give the equation (2) above ? 

9.3 Put the equations with matrix notation 

(i) 3x 2 4- 2 xy +y 2 — 1 

(ii) 4x 2 - 3xy + 2y 2 = 1. 


9.2 Reduction to Diagonal Form 

The general homogeneous quadratic form ax 2 4- 2 hxy + by 2 = c can be 

represented in matrix form asx' ^ x=c where the matrix ^ ^ 

has been chosen here as a symmetric matrix - the most useful form, as 
will be seen from later work. 

This quadratic form represents a conic section, which will be a pair 
of straight lines, a circle, an ellipse or a hyperbola depending on the 
values of a , b and h. The centre of the conic is the origin but its axes 
do not necessarily lie along the coordinate axes. It is possible to find the 
direction of the axes of the conic and its form, by finding the eigenvector 
of the corresponding matrix. 

Consider the quadratic form 6x 2 + 4 xy + 3 y 2 = 14. When written in 
matrix form x'Ax = 14, with A a symmetric matrix, we have A = ^ ^ • 

The eigenvalues of A are given by 

(6 - A)(3 - A) - 4 = 0 
A 2 - 9A +14 = 0 


i.e. 

i.e. 

and 


(A - 7)(A - 2) = 0 

A = 2 or 7. 

When A = 2, from (6-A)# + 2jy = 0 the eigenvector is given by ^ }^J. 

When A = 7, the eigenvector is given by ^ ^ (perpendicular to the first 
eigenvector, as A is symmetric). 

Considering the eigenvector with positive gradient, ; if we make 

this an axis, then the matrix which will rotate the axes to this position is 

given by R = ( C . 0S a Sm a j where tan a = A, as for the eigenvector 
8 J \ - sin a cos a J 8 

;#= 2 , jy = l;soR=-^^ ^ 2 ^- What will happen to the quad¬ 

ratic form when we refer it to axes taken with respect to the eigenvector 
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of the matrix? If the new axes are labelled x 1 y 1 then the rotation matrix 
has the effect x x = Rx where (x v y x ) are the coordinates of the point 
(x y y) when referred to the new axes. 

From x 2 = Rx then x = R -1 x r 

Hence the original quadratic form x'Ax = 14 becomes 


(R _1 x 1 )'A(R _1 x 1 ) = 14. 

Hence x^R-^'AR^x, = 14 as (R-^x,)' =x^(R- 1 )' or x^RAR^ = 14 as 
R _1 = R' and (R')' = R. The quadratic form considered thus becomes 



(-i 2)(2 3)75(1 2) x i- 14 


The matrices reduce to 



- a matrix which is said to be of diagonal 


form, that is with entries only in the leading diagonal, other values being 
zero. The values of the entries in the leading diagonal of this matrix are 
the eigenvalues of the original matrix. 


The quadratic form has been reduced to 


4 °) xi=14 

^ 4 2 ) (;;H 

7xf + 2yl = U, 


an ellipse with semi major axis of J2 and semi minor axis of J7. 


- 



Fig. 9.1 
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9.3 Worked investigations 

(i) Interpret the form of the quadratic equation 4# 2 -12 xy 4- 9y 2 = 13. 
The quadratic form in matrix notation is x' ^ ^ ^ x = 13 where 

the matrix is symmetric. The eigenvalues of ^ ^ ^ are given by 

(4-A)(9-A)-36 = 0 


A 2 - 13A=0 

A = 0 or 13. 


The eigenvector when A = 0, from (4 - A)* - 6y = 0 is an< 3 the eigen¬ 
vector when A = 13 is ^ ^ . 

Considering axes rotated so that the eigenvector becomes the new 


Xj axis; the rotation matrix is given by 


s/13 


iUJ d - 


as tana 


and thus the quadratic form when referred to x 1 y 1 axes becomes 

, i 


V13 


u du D4,a I)— 


The matrices reduce to ^ 13 ) ” aga * n f° rm anc l w hh 

diagonal entries equal to the eigenvalues. 

The quadratic form now is 

x (0 is ) 11 - 13 

<*■ *>(S 1 X:)- 13 ' 

13 ^ = 13 , 

;Vl=±1 

which represents two straight lines. 

(The original equation was Ax 2 - 12xy + 9y 2 = 13 giving 
2x-Zy = 4-s/13 or 2x -3y = -s/13; two straight lines.) 


or 


or 


and 
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Note also that the vanishing of the determinant of the matrix, which 
implies an eigenvector is zero, is also the condition for the second degree 
terms to form a perfect square. 

fii) Interpret the form of the quadratic equation 
x 2 - 8xy + ly 2 = 9. 


In matrix notation this becomes 


■•(-I -})*-»■ 


The eigenvalues are given by 

(1 -A)(7 -A) -16=0 
7-8A+A a -16=0 
A 2 -8A-9=0 
(A-9)(A + 1)=0, 

A =9 or 


therefore 


• 1 . 


The eigenvector, when A = 9, from (1 - A)* - 4y = 0 is ^ ^ an d when 
A *= -1 the eigenvector is . Taking the eigenvector , of gradient £ 
the rotation matrix is ^ • 
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The quadratic form when referred to the eigenvector as new x 1 axis is 
given by 




or 


(-1 2) (-4 7)75(1 2) Xi ~ 9 

X2 ( o 9) x 2 = 9 ’ 
-x\ + 9y\ = 9, 


- +yf = 1, a hyperbola. 


Examples 

Investigate and sketch the following quadratic forms as conics. 

9.4 3x 2 -2xy + 3y 2 = 8. 

9.5 x 2 + 4xy-2y 2 = 6. 

9.6 4x 2 -4xy+y 2 = 20. 

9.7 4x 2 + 4y 2 = 9. 

An equation of the form ax 2 + 2 hxy + by 2 = c cannot be a parabola as a 
change in sign of x and y leaves the equation unaltered, and hence it 
represents a curve symmetrical about the origin. 


9.4 Formulation and proof of general results 

In the examples which have been worked through, several key points 
always appear. These need to be looked at generally and proofs of them 
established. 
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The quadratic form used viz: 6x 2 + 4 xy + 3 y 2 = 14 was put into matrix 
notation asx'^ ^ x = 14 with the symmetric matrix ^ 3 ^ • The 
eigenvalues of this matrix were found to be 2 and 7 and the correspond¬ 
ing eigenvectors ^ _ 2 ^ ^ * When the quadratic form was referred 

to the eigenvector as axis, its resultant form became 2) ~ 

matrix of diagonal form with entries in the leading diagonal of the eigen - 
values , with the eigenvalue 7 (corresponding to the eigenvector used) in 
the top left hand corner. Will these results be true generally? 

Consider a symmetric matrix A which has eigenvalues of X x and A 2 
and which will have eigenvectors which are orthogonal. If we choose 
one of the eigenvectors to be an axis and rotate the axes according to the 

rotation of axes matrix R = ( C . 0S a Sln a ] then the matrix A trans- 

\-sina cos a ) 

forms into RAR' (see page 57). 

The eigenvalues of RAR' have been found to be and A 2 and RAR' 


is of the diagonal form ^ ^ ^ . 


The eigenvalues of RAR' are given by the roots of the equation 
det (RAR' - AI) = 0. 

But RAR' - AI = RAR' - RAIR' since RI=IR and RR =RR 1 = 1 
= R(AR' - AIR') by the distributive law 
= R(A - AI)R' by the distributive law 

therefore 

det (RAR' - AI) = det R det (A - AI) det R' (page 30) 

= det (A - AI) since det R = det R' = 1 

that is A and RAR' have the same characteristic equations. These 
equations must therefore have the same roots, and so have the same 
eigenvalues. 

To prove RAR' = ^ ^ ^, let one eigenvector be given by the unit 
vector c = and the other eigenvector by the unit vector d = • 

Then ( 1 ) may be written as ( C . 0S a ) for some a and (j | will be 

W \ sin a / W 

represented by ^ ^ since c and d are orthogonal as A is symmetric 
(page 53). 



120 


MATRICES AND THEIR APPLICATIONS 


axis 


Then c'c = (c x c 2 ) (^j = c x + c\ = cos 2 a + sin 2 a = 1 and similarly 

d'd=if + (/l = l.( 1 ) 

c'd = (c x c 2 ) + c 2^2 = cos a . - sin a + sin a . cos a = 0 = d'c (2) 

R is the rotation of axes matrix needed to make the eigenvector c the x 

is; hence R = f C .° Sa Sm a W C J C A. 

’ \ - sm a cos a / d 2 ) 

'=fa c Aa( Ci d A 

\d x d 2 J \c 2 d 2 ) 

-(X X ) <** 

*a\ (^i c A 2 d) since Ac=AjCand Ad=A 2 d, c and d 

</j d 2 ) being eigenvectors, 

D 


SoRAR 


_ f C l r 2^\ (\ C 1 A 2 J X 

\^1 d 2 ) \A^2 A 2^2 

_ Aifci + c 2 ) A^^ + c 2 d^\ 

\X^(c 1 d 1 + c 2 d^) X 2 (df + d 2 ) ) 

“(o 1 a 2 ) usin g (!) and (2). 


9.5 Eigenvectors and Principal Axes 

We have seen that when the axes of coordinates are rotated into the 
positions of the eigenvectors of the matrix A, the equation x'Ax = & 
reduces to the form A^f+A^y 2 ^ which is the equation of a conic 
(possibly degenerate) with its principal axes along the axes of coordinates. 
Thus we may say that the directions of the eigenvectors of A are the 
directions of the principal axes of the conic x'Ax = k. 


9.6 The general quadratic form - translation of axes 

We have, at present, only considered quadratic forms containing no 
first degree terms. The most general form is given by 

f(x,y) = ax 2 + 2hxy+by 2 + 2gx + 2fy+c . . . (3). 

We have seen how the part ax 2 + 2 hxy + by 2 can be reduced to diagonal 
form by a rotation of axes. What effect have the terms 2 gx and 2fy on 
our previous analysis? 
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Consider the form 

x 2 +4xy-2y 2 -10x + 4y+20=0 . . . (4). 

We can use matrices to describe this form also. Check that 


/ 

' 1 

2 

“ 5 \ 

fX\ 

(x y 1)1 

2 

-2 

2 ) 

U Uo 

\ 

is the same as equation (4). 

<-5 

2 

20 / 

Vi/ 


( 1 2 ~ 5 \ . 

Note that the matrix A = ( 2 -2 2 1 is again chosen to be 

V-5 2 20/ 

symmetric. 

Show that the general quadratic form (3) can be written as 

(a h g\/x\ 

(* y 1)(A b /)(*). 

Make up other.forms of the matrix A which will be equivalent to (4). 



We shall see that the terms 2gx, 2fy can be removed by a translation 
of the axes, and that the resulting expression can then be reduced to 
diagonal form as in the preceding section. 

The origin is to be translated from O to O '(r y s). If the new axes are 
x 3 y v we have 

x 1 =x-r 
y x =y-s. 

These equations in matrix form are 
*i\ /I 0 

*H° 1 

1 / \o 0 

or Xj =Tx where 

/I 0 
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If we use this translation what will be the new form of the expression 
x'Ax? 

By using x = T^Xj, x'Ax becomes 

(T-i Xl )'A(T-i X ) =xJ[(T-i)'AT-i)] Xl = X ;Bx x say. 

/I 0 r\ /I 0 0\ 

T-!= 0 1 * and (T-i)'=(o 1 0 . 


0 0 1 


r s 1 


( a h g\ 

For the matrix A = I h b / ) (T - 1 )'AT -1 becomes 

V / c) 

A 0 0 \/a h g\ A 0 r\ 

0 1 0 (A b f 0 1 5 

v s 1 / f c) Vo 0 1 / 

( a h ar + hs +g \ 

h b kr + bs +f j 

ar + hs +g hr + bs+f ar 2 + 2 hrs + bs 2 + 2 gr + 2fs + a 

( a h ar + hs +g\ 

h b hr + bs+f J = B. 

ar + hs+g hr + bs+f f(r y s) / 

The symmetric form of this matrix is a consequence of taking A to be 
symmetric. 

We wish to choose r and s so that the form of B becomes 
px\+2kx 1 y l + qy\+c v 
This expression in matrix form is 

(P k 0\/*A 

(*i yi ? °)(>'i)- 

Vo o cj \i / 

Hence we chose/) = a y q = b,k=h,c 1 =f(r, s) and 

ar + hs+g = 0 . . . . (5) 

hr + bs+f =0 . . . . ( 6 ) 

/ a h g\ 

(Note the first row of A = ( A b /lisa h g and the second row 

V / cJ 

h b /.) 

The need for A to be symmetric is that the four zeros in B are all 
produced. 
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The equations (5) and ( 6 ) allow us to find the values of r and s , so that 
the quadratic form is reduced to ax f + 2hx 1 y 1 + by\ +/(r, s). By methods 
described in the first section, this form may now be reduced to 
mx\ + ny\ 4 - c 2 by a suitable rotation of axes and hence the locus of 
f(x, y) = 0 determined. 


9.7 Worked investigations 

(1) x 2 + 4 xy -2 y 2 -10x + 4y+ 20=0 =f(x,y) can be put in matrix 
form as 



Let A=( 2 -2 2 1. Then by translating the axes so that 

V-5 2 20/ 

(r, s) becomes the new origin we must have (for the quadratic form to 
reduce to ax f + 2hx 1 y 1 + by\ + c x = 0) that 

lr+ 2y -5 =0 

and 

2r - 2s + 2 = 0 

from which r = 1, s = 2. 

Then/(1, 2) = 1 +4(2) -2(4) - 10(1) +4(2) +20 = 19. Hence by trans¬ 
lating the axes to the point (1,2) the equation becomes 

/I 2 0\ /*A 

(*i y x 1) 2 -2 0 U =0 

Vo 0 19 /Vl/ 

or 

<- J)(5) +19 - 0 - 

The form (*, J.lQ (*|) can now be reduced to mx\ + ny\ by a 

suitable rotation of axes. 


that is 


Let C = ^ 2 ) * eigenvalues of C are given by 

(1 -A)(- 2- A)- 4 = 0 or A 2 +A -6 = 0=(A + 3)(A-2) 

A = - 3 or 2. 

From (1 -X)x + 2y = 0> when A= - 3 the eigenvector is anc * 

when A =2, the eigenvector is (as expected, perpendicular to the 
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first eigenvector, as C is symmetric). The rotation matrix to rotate the 

axes through tan -1 1 to make the eigenvector the axis x 2 is 


U-l i)- 


Hence the matrix C under rotation become RCR', or 

l ( 2 1\/1 2 \ J (2 - 1\_ a /4 2\/2 - 1 \ 

M" 1 2/ \2 2j-*[3 -6/ \1 2) 

_i/10 0W2 0\ 

5 V 0 -15) V° - 3 / 

(Note the eigenvalues in this diagonal matrix and the eigenvalue 2 


corresponding to the eigenvector ( 1 ) taken.) 


The equation (aq y x ) 


U 


+ 19=0 when referred to axes 


of x 2 , y 2 inclined at tan -1 (|) to the x l9 y ± axes, becomes 

('• +o -S)(?J + I9 -°- 

2 * 1 - 3^1 + 19=0 

or 2*1. 3^5, ,_. 


- t — 

19 19 


= 1 - a hyperbola. 



Fig. 9.5 
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(2) Determine the locus of the conic 

4x 2 + 2 xy + 4y 2 - 26x - 14 y + 31=0. 

/ 4 1 

f(x,y)=4x 2 + 2xy + 4-y z -26x-14y + 31=x'l 1 4 

V —13 -7 

( 4 1 ~ U \ 

LetA = ( 1 4 -7 . 

V —13 -7 31/ 



x = 0 . 


Consider a translation of axes so that the point (r, $) becomes the new 
origin. The values of r and s are given by 


4r + s - 13 = 0/ 
r + 4s- 7 = 0/ 


from which r = 3 and $ = 1. 


Then/(r, s)=f(3, 1) = -15. 


On referring the conic to the origin (3, 1) with axes its form 
becomes 

A 1 (k 
xi(l 4 0 Jx x =0 

\0 0 -15/ 

or 

X ’(l 4 ^ x-15==()i 

Let B = ^ 4 * ). The characteristic equation of B is given by 

( 4 i A .V) - 

that is 

(4 -A ) 2 -1=0 
A 2 - 8 A + 15 =0, 

A = 3 or 5. 


When A = 3 the eigenvector is ^ and when x = 5 the eigenvector 
is ^. Rotating the axes to the eigenvector , through 45° the form 

of the conic becomes *2 3 ^ x 2 -15 = 0 . 

(Note that the eigenvalues have been put in the diagonal matrix with¬ 
out working out the matrices with the rotation matrix through 45° and 

that since the eigenvector ^ was chosen corresponding to eigenvalue 5, 

this eigenvalue takes the top left hand corner.) 
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is 


* 2 (o a)* 2 15-0 

5*1 + 3?!-15=0, 

*1 y\ < 

^= + — = 1 , an ellipse. 



Fig. 9.6 


(3) Determine the locus of the conic x 2 + 4 xy + 4y 2 - 28x - 6y - 4 = 0. 

/ 1 2 - 14 \ 

/(*,?)=* 2 +4*y+4? 2 -28*-6?-4=x'( 2 4 -3)x = 0. 

V-14 -3 -4/ 


/ 1 2 " 14 \ 

Let A=( 2 4 -3). 

V-14 -3 -4/ 


Consider a translation of axes so that the point (r, 5 ) becomes the new 
origin. The values of r and s are given by 

r+ 2$-14 = 0 
2r + 4*:- 3=0 

which are inconsistent equations and hence no values of r and s. 

We are thus unable to translate the axes to a suitable origin. 

Let B = Q ^ . The eigenvalues are given by (1 -A)(4 - A) -4 = 0, 


- 5A = 0; A = 0 or 5. When A = 5 the eigenvector is 


Rotate the axes through tan -1 2 . The matrix B will become 
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but the values of the other terms of A will also be affected by the rota¬ 
tion. In this particular case when there are no (r, s) values, we need to 
see more closely the effect of the rotation. The rotation matrix is 

^ _ 2 j ) fr° m axes x y t0 axes say; so 

(;iBU ?)(;) 


and hence 


so that 


CHG D(;:) 






In/(#, y) the terms - 28# - become 


(- 28# x + 56^ -12#! - 6y x ) 


- 75 (-*0 Xl+ 50 yi) 
= - 8n/5#! + IOn/Sv!). 



Fig. 9.7 
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The conic f(x , y) = 0 then becomes, on rotation of axes through tan -1 (2); 


* (o J) X] - 8^5*, +10 J5y x -4 = 0 


or 


+ 10^5)^ -4 =0 

- which is a parabola, since the expression contains no y\ term. 


— ■ 


V5 3 


—t - zJsyn 


(*i J) a =t-2v/5j 1 +^=4-2V5y 1 . 


-4 


Let z 1 = x i and z 2 = lj5y 1 - 4 then z\ = - z 2 , 


9.8 The geometrical significance of the invariants in the matrices 
used in the reduction of quadratic forms. 

If x'Ax, where A = ^f is a reduced quadratic form, then the 

geometrical characteristic of the locus x'Ax + c = 0 are given by the 
following invariants of the matrix. Let the eigenvalues of A be X 1 and 
A 2 then, for X x and A 2 , 


from which 


(a-X h \ Q 

l h b-\) 

(<a-A)(6-A)-A 2 = 0, 
A 2 -(a + 6)A + «&-A 2 = 0. 


Hence 

and 


Aj + A 2 = cl + b 
AjA 2 ~ab-h 2 = det A. 

On using a rotation matrix to make the eigenvectors the axes the 
diagonal form of A is ^ ^ ^ , that is the quadratic form reduces to 

x i(o A 2 ) Xl+C = 0or 


AjXj + A 2< yf c = 0 . 


. (7). 
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Now if in (7) c — 0, the quadratic form is a line pair going through the 
origin of the axes a?j y v If A 2 +A 2 = 0 the lines are the orthogonal pair 

±*i- 

In other cases from the form of (7) we have 

(i) an ellipse if Aj and A 2 have both the same sign, that is A X A 2 > 0, 
that is ab-h 2 > 0 

(ii) a hyperbola if A x and A 2 are of different signs, that is A 1 A 2 < 0, 
that is ab - h 2 < 0 (rectangular when A 1 -f- A 2 = 0) 

(iii) a parabola if one of A x or A 2 = 0, that is ab - h 2 = 0 

(iv) acircleifA 1 =A 2 . 

(Note that it is possible for (7) to give imaginary conics.) 

Examples 

Determine the locus of the following conics and sketch them 

9.8 5x 2 + 4 xy + 5y 2 + 6x - 6y - 15 = 0. 

9.9 4a? 2 + 4 y 2 - 16a? - ly +4 = 0. 

9.10 2a: 2 +4xy -y 2 + 4x - 8y -16 = 0. 

9.11 a? 2 - 2xy +y 2 -2x -2y = 0. 
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ANSWERS TO EXERCISES 


Chapter 1 

1.1 ( 7sV 35 V 40^ (145,120, SO, 315), (-10, 30, - 5,15). 
\170/ V 30/ M05/ 


1.2 


e 

a 

f 

r 


e 

a 

b 

*3 


m 

25 

30 

0 

55 

m 3 

25 

30 

1 

54 

s 2 = 

w 

35 

18 

0 

53 

S 3 = w 3 

35 

18 

2 

51 


c 

60 

48 

0 

108 

c 3 

60 

48 

3 

105 

1.3 



e 

a 

f 


e 

a 

b 

*3 

F- 

-S = 

m 

20 

20 

15 

CO 

1 

Xfl 

CO 

II 

H 

CO 

20 

20 

1 

39 



w 

20 

2 

20 

W 3 

20 

2 

1 

21 







C 3 

40 

22 

2 

60 


e a f 

4S= m 100 120 0 

w 140 72 0. 


1.4 English students’ subscription =45.3 +55.2 = 245=(45, 55).(3,2) 
Maths students’ subscription = 50.3 + 20.2 = 190 = (50,20) .(3,2 
All students’ subscription =93.3 + 72.2 = 423 

(this is the true total). 

1.5 -6,0,1 

1.6 ab_(j "),ba,(j 

“-(J ?)’ da -G 1 ).bc-cb=b.bd.(^ J). 
DB -(l "),CD-DC.D ; A>-(” “),B*-(J 2), 

Hi ?)• 
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MATRICES AND THEIR APPLICATIONS 


1.7 EF does not exist (for matrix multiplication to be possible the 
number of columns of the first matrix must be equal to the number of 
rows in the second matrix). 

**~{l 6 "is) EG-GE-E, EH.( j 1 -j) 


/I 

1 

-!\ 

HE= 2 

3 

5 1, FG=F, GF does not exist, 

Vo - 

2 

1 / 

™-(5 3 

-1 

3 

,), HF does not exist; GH=HG=H. 

/ 7 

-1 

12 \ 

E 2 = 3 

6 

3 ). F 2 does not exist. 

V -2 

-4 

3/ 


C and G leave matrices unchanged; these are called unit matrices. 

D and H when pre-multiplying a matrix interchange its first and 
second rows, when post-multiplying a matrix interchange its first and 
second columns. 


1.8 


/ 

t = Mv where v = 

\ 


1 

1 

1 

1 

1 

1 


gives the t column 


/ 7?\ 
131 
128 
207 
\230/ 


= 713; 


t s v = 924 and as 1*30 = fff the value 1*30 comes from using mvt c = t s . 


/i\ u=a 

1.9 Let N be the 6x5 matrix, v = 1 I, 

l;/ '- p 


iiiii), 
2 3 4 5 6). 


Nv gives the number of families with 1 child, 2 children; the total 
column. 

uN gives the number of families in each group. 
pN gives the number of children in each group. 

If m represents the row giving the mean number of children in each 
group, then mv(uN)=(pN). 
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Chapter 2 

2.1 (a) fa 0\ (b) (k 0\ (c) / -1 0\ 

VO l) [0 k) \ 0 l)‘ 

2.2 (a) reflection in x axis 

(b) linear stretch, factor 4 in # axis direction 

(c) enlargement, factor 3 

(d) rotation of the plane through 45° 

(e) rotation of the axes through 45°. 



2.5 (i) A shear parallel to the x axis in the negative direction, factor 1. 

(ii) A shear parallel to the axis of factor 
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MATRICES AND THEIR APPLICATIONS 


(c) 



2.7 Identity, reflection in y axis, reflection in x axis, rotation through 
180°. 

Reflection in x=y, rotation through 90°, rotation through 270°, 
reflection in x = -y. 


2.8 



The matrix ^ ^^ transforms the plane such that P(x , y) is mapped 

onto P*(x* y y*) where (^y*J = (^j • Then x* =2x and y* =y. 

Hence if (x y y) are related by the equation x 2 +y 2 = 1, then (x* y y*) are 
/x*\ 2 

related by (— 1 + O '*) 2 = 1, which is an ellipse. 

2.9 (a) Shear factor 2 parallel toy axis. 

(b) Projection parallel to x axis onto linejy = 3x . 

2.10 /\ A 

V2\i -1/ 





ANSWERS TO EXERCISES 


135 


2.11 M - reflection in y axis; Q - rotation through 90°. 

QM = ^ j - reflection inx= -y. 

MQ = ^ ^ - reflection in x =y. 

t\ =A 

2.12 R = l i I > find th e new positions of the vertices of the 

\T 2 / 

unit square and plot them. 

S = ^ ^ , use the coordinates of the rotated square and find their 

new positions due to the shear. 

1 3 -V3 X 

2 


SR = | 


I, find the new positions of the vertices of the 


2 2 
J3 3J3 1 

2 2 + 2 

original unit square and plot them. The final figures will be the same. 
1 - 1 \ / 1 -1 


- 1 ) f 

V2 J2‘ 

i)d -?)f 

J2 / \J2 

( c ) /0 1\/1 0\/2 0 W 0 6 \ 

\i o/\o 3/ \o 2)-\2 o;- 




2.13 

(a) 




B 




-*~X 


rotation through 270° 
(note position of vertices) 



MATRICES AND THEIR APPLICATIONS 


%.:)(! ;h-? o)- 


^ ^ /cos a -sina\ ^ _/cos/3 -sinj8\ 
a \sina cos a/ P~\sin P cos pj 

P _ /cos (a + j8) -sin(a+j8)\ 

a \sin (a + /?) cos (a + /?)/ 

P P _ /cos p - sin j8\ /cos a - sin a\ 

^ a \sincos p) \sin a cos a ) 


2.15 A-(; -J) 

i) 

M -(°l J) 
MN -(-? J) 



( q jx results the same: trans- 
_ j q J formation of unit square 
' would correspond. 



/q _ jx results different, transfor- 
NM = f ^ J mation of unit square would 
' ' not correspond. 


AB = MN: a reflection in two mirrors is equivalent to a rotation 
through twice the angle between the mirrors. (MN implies reflection 
in y axis first followed by reflection in x =y.) 


2.16 AB = BA and A(BC) = (AB)C. 


2.17 (a) (2 4\ (2 

V> 2) “VO 

S)G *) 

a shear followed by an en¬ 
largement (or vice versa) 

(b) (? J) -(? 

i)Co ?) 

a shear followed by a reflec¬ 
tion in x =y (not reversible) 

^kTo 

1 

1 ! 

-l\/2 O' 
0 A° 3, 

a double stretch followed by 
) a rotation through 90° 
(reversible) 

(-! ?)-(- 

1 0 \/l 0 > 
0 1/ \3 lj 

i a shear followed by reflection 
' in y axis (not reversible). 
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2.18 A-(J J) B-(J _J) 

AB = rotation through 90° 

BA = rotation through 270° 

AB 2 =A 

(BA) 2 = rotation through 180° 

ABA = reflection in y axis 
BAB = reflection in x — -y 

ABABA=A(BA) 2 = ( _ J ” J j = BAB 

BABABA=B( ABABA) = B(BAB) = AB. 


2.19 


2.20 




I 

R 

+ 

0 

1 






I 

1 

h 

~0~ 

0 

”1 






R 

R 

I 

1 

1 

0 





V 

-( 

-1 

0 

?) 

H = 

-G 

1 

3 

0> 

-b 

) 

R 

-I 

I 

-( 

1 0 
0 1 

) 









I 

R 

V 

H 


X 

1 

3 

5 

7 

T~ 

I 

R 

V 

H 


1 

1 

3 

5 

7 

R 

R 

I 

H 

V 


3 

3 

1 

7 

5 

V 

V 

H 

I 

R 


5 

5 

7 

1 

3 

H 

H 

V 

R 

I 


7 

7 

5 

3 

1 


(Mod 8) 


Chapter 3 


U W ( 7 1) (ii) ( _ 7 “1) (iU) ( 5 " 3 ). 


Interchange the entries in the major diagonal and reverse the sign of 
the entries in the minor diagonal. 


3.2 (i) 

-'•( 

1 A) 

(ii) -2, ^ 

-2 

3 

2 

(in) 

-5, i 

Cl f) 

(iv) 13, * 

(-5 

(V) 

/ cos a 

-sina\ 




\sina 

cos a) * 
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3.3 


MATRICES AND THEIR APPLICATIONS 




The matrix A is a non-singular matrix which has the geometric 
property of 1 to 1 correspondence between points and their images. The 
matrix B is singular and gives a many-to-1 correspondence between 
points and their images. For B any point on the line x + 2y = 3 is 
mapped onto (3, 6). 

3.5 /cos a - sin a h cos a 

I sin a cos a h cos a 

V 0 0 


- k sin a\ 
+ k sin a I 
1 / 



3.8 (i) /I 0 0\ (ii) (k 0 0\ (iii) /I a 0\ 

0 10) (0 A 0) [b 1 0). 

\0 0 -1/ Vo 0 k) Vo 0 1/ 
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3.9 (i) reflection in the x - z axes plane 

(ii) a linear stretch of factor a in the x axis direction 

(iii) identity 

(iv) shear factor 1 parallel to z axis. 


Chapter 4 

4.1 (a) x=2,y = 3; (b) x = 3,y = 2; (c) x= -},y = If. 


4.2 

(i) 

x = 2,y = l; 

(ii) x 

=2,y=- 

i; 

(iii) x = 

hy 


Civ) 

£ 

II 

O 

II 

H 

(v) * 

= l,y = 3. 




4.4 

(i) 

( 4 3 

“ 3 \ 

(») i/ 

-5 

7 




-4 -2 

3 ) 


-1 

-1 

1 



V-i -1 

1/ 

\ 

-3 

-5 

1 / 


(iii) 

|/-5 9 

-3> 

i (iv) / 

5 

-7 

-2 



3 -3 

-3 

/ 

' 

7 

-10 

-3 



\-7 9 

-l) 


-is 

22 

6 





\ 

-3 

4 

1 


Chapter 5 

5.3 (a) (i) 1, (ii) 2, (iii) 0, (iv) -3, (v) 1. 
(b) (i) and (v). 

5.4 (i) 1,6; x-4y = 0, x+y = 0; 

(ii) 5,10; x+2y = 0,2x-y = 0; ( _?). (J) • 

5.5 2, 1; ± i . 


5.6 A 2 -X(a-\-d) +ad-bc—0 f 

Sum of eigenvalues is sum of major diagonal entries, product of 
eigenvalues is determinant value. 


5.7 See 5.5(ii). 


5.8 



has repeated eigenvalues of 2. 
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5.9 Eigenvalues are major diagonal entries. 

5.10 k and k . 


5.11 a±b. 

5.12 7,2; 7,2; -5,8; -5,8. 

5.13 x + 3y =0, x-y = 0; x axis, x=y; x + 2y = 0, complex. 

5.14 (a) real distinct eigenvectors 

(b) complex eigenvectors 

(c) any line through the origin is an eigenvector 

(d) one eigenvector 

(e) x zndy axes as eigenvectors. 

5,5 (-;)■(:!)■(:!)• 

(ii> -I. 2 .1; (o). (]). 0- 



5.17 Eigenvalues 3 or - 7, eigenvector lines x - 2y = 0, 2x + b = 0. The 
eigenvectors are orthogonal. 

The characteristic equation for M = ^ ^ is 

A 2 -A(a + &)+a&-c 2 = 0. 

For the equation to have real and distinct eigenvalues 


(a + &) 2 - A(ab - c 2 ) >0. 

The left hand side is (tf-&) 2 + 4c 2 which is greater than zero for all 
a , b , c. If u and v are the eigenvectors with corresponding eigenvalues 
Aj and A 2 , then Mu = A x u and Mv = A 2 v. 

From Mu =A x u, then (Mu)' =X 1 u' or u'M' =A 1 u', since (Mu)' =u'M' 
(see 5.16). As M = M', u'M=A x u'. Postmultiplying by v, u'My=A 1 u'v 
and as Mv=A 2 v, u'A 2 v = X{a'y. Thus A 2 u'v - ^u'v = 0 or (A 2 - A x )u'v = 0 
and since A 1 ^A (eigenvalues are real and distinct), u'v = 0, showing the 
eigenvectors u and v are orthogonal. 
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/ 4 3 3\ 

5.18 M- 1 - -4 -2 3). 

V-l -1 1/ 


7.1 



Chapter 7 

a & £ d e 

a 1 0 1 1 0 1 \ 

*110100 
c 1 1 0 1 1 . 

<f\0 0 1 0 1 

* \l 0 1 1 0/ 



The matrices are symmetrical since if <2 is joined to b is joined to a . 
The sum of the individual rows (or columns) give the number of 
connections to (or from) a node. 

When the matrix is not symmetrical the sum of entries in a row give 
the number of connections from that node to other nodes and the 
sum of the entries in the columns give the number of connections to that 
node. 



7.4 /I 0 1\ /0 1 1\ /I 1 2\ 

1 0 1+0 1 01 = 1 1 1 ). 

Vo 2 0/ Vo 0 0/ Vo 2 0/ 
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7.5 (a) the number of links to a station 
(b) each link connects two stations. 

7.6 Rows and columns containing two l’s-rows b , c and columns 3,4; 
also rows a , b> c and columns 1, 2, 3. 

7.7 (a) 1 and 2 from the a row. 

(b) Links 2, 3 and 4 cut, but not 1. The links to cut are found by 
adding modulo 2 (so 1 -b 1 =0), the entries in the columns of 
the rows a and b , and cutting that link when the result is 1. 

abed abed 

1/1 1 0 0\ a /2 1 1 0\ 

2/1 0 1 0\ b 132 01 

7.8 A' = 3 0 1 1 0 AA = c\l 2 3 0/ 

410 1 1 0/ d\ 0 0 0 1/ 

5\0 0 1 1/ 

The entries in AA' give the number of links between different 
stations; the leading diagonal giving the number of links at that station. 

a ft y a |3 y 

1 /I 0 1\ /I 0 1\ 

2/1 0 ll b 1 1 l\ 

7.9 P= 3 1 1 0 Q = c\l 1 1 1 

4 0 11/ d\0 0 1/ 

5 \0 0 1/ 

1 2 3 4 5 

1/221 1 1\ The entries are the number of 

2/221 1 1 regions common to the links; the 

PP' =311210 leading diagonal giving the number 

411 1 1 2 1 of regions with each link. 

5\1 1 0 1 1/ 

a P 7 

a /3 1 2\ The entries are the number of links com- 

P'P= 1 2 1 j mon to the regions; the leading diagonal 

y^2 1 4/ giving the number of links with each 

region. 

abed 

a 12 2 2 1\ The entries give the number of 

QQ' = b 12 3 3 11 regions common to the stations; the 

c 1 2 3 3 1 / leading diagonal giving the number of 

d \l 1 1 l/ regions with each station. 
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(ii) not possible, since the column of 
each link must have two, and only 
two l’s indicating the two connected 
stations. 



a 

b 

c 

d 

e 


a 

b 

C 

d 

e 

a 

(° 

0 

0 

1 

°\ 

a 

1 ° 

0 

0 

0 

i\ 

b 

'o 

0 

1 

1 

o' 

b 

[i 

0 

0 

0 

i 

c 

1 

0 

0 

0 

0 

M 2 = c 

0 

0 

0 

1 

0 

d 

1 ° 

0 

0 

0 

ll 

d 

[1 

1 

0 

0 

0 

e 


1 

0 

0 

0 ) 

e 

\0 

0 

1 

2 

0/ 


Powers (M + M 2 ) 
a 2 
b 4 
c 2 
d 3 
e 5 

e has the most 1 and 2 stage connections with the other stations. 



(ii) A 2 = 


a b c d e 

a 1 0 0 1 1 1\ 

b 1 0 2 0 l\ 

c 0 10 12 

d 12 10 1/ 

e\l 1 0 0 0/ 


(iii) Yes. (iv) A + A 2 = 



a 

b 

c 

d 

e 

Powers are 

a 

10 

1 

1 

1 

2 \ 

a 

5 

b 

1 

0 

2 

1 

2] 

b 

6 

c 

1 

2 

0 

1 

2 

c 

6 

d 

2 

2 

2 

0 

2 J 

d 

8 

e 

\l 

1 

1 

0 

0/ 

e 

3 


d can communicate with all others (twice) in one or two stages. Only 
e cannot communicate (with d) in one or two stages. 
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MATRICES AND THEIR APPLICATIONS 



M = 


A B 
A /0 1 

ft / ( 


C 

D 

E 

F 


0 0 
1 0 
0 0 
0 
0 


V! 


C 

1 

0 

0 

1 

1 

0 


D E F 
0 0 1 
0 0 0 
1 1 0 
0 1 0 
0 0 1 
0 1 0 


A B C D E F Powers M+M 2 are 

A / 2 0 0 1 2 0\ A 8 

B/0 000 0 0] B 0 

M 2 = Cl 0 1 3 0 1 2 C 10 

D 10 1111 D 7 

E \ 2 0 0 1 2 0 E 7 

F \0 1 2 0 0 2 / FI 


Mrs. C will spread a rumour quickest throughout the group. 


7.14 


A = 


R = 


A 2 — 


a 

b 


c 

d 


e 

f 


a 

b 


c 

d 


f 


a 

b 


c 

d 


f 


( 


a 

b 

c 

d 

e 

/ 

0 

1 

0 

1 

0 

r 

0 

0 

1 

1 

0 

i 

0 

1 

0 

1 

1 

l 

1 

1 

1 

0 

0 

i 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0, 

a 

b 

c 


d 


0 

0 

-1 


0 


0 

0 

0 


0 


0 

0 

0 


0 


0 

0 

0 


0 

- 

1 

0 

0 


-l 


0 

0 

0 


0 

a 

b 

c 

d 

e 

/ 

2 

2 

3 

2 

0 

2 

2 

3 

2 

2 

1 

2 

2 

2 

3 

2 

0 

2 

1 

3 

2 

4 

1 

3 

0 

0 

0 

0 

0 

0 

1 

3 

2 

3 

1 

4 


o 

-l 

0 

-1 

0 

0 
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A + A 2 + R = 


a 

b 


c 

d 


e 

f 


a 

b 

c 

d 

e 

/ 


2 

2 

2 

3 

0 

3 \ 

12 

2 

3 

3 

3 

0 

3 i 

14 

2 

3 

3 

3 

1 

3 

15 

2 

4 

3 

4 

0 

4 

17 

1 

0 

0 

-1 

0 

0 ) 

-2 

2 

4 

3 

4 

1 

4/ 

18 

9 

16 

14 

16 

2 

17 



Social index a = 9+2(12) =33 
b = 16 +2(14) =44 
c = 14 + 2(15) =44 
d= 16 +2(17) =50 
e = 2 + 2(-2)= -2 

/= 17+ 2(18) =55. 


Chapter 8 

8.1 3,000-1st year; 6,916 (approx.) 2nd year; =5,666; 3,000; .... 

8.2 (a) f of the species survive their first birthday; of these J survive 
their second birthday and live into their third year and none of the 
species live longer than 3 years. In the third year 6 of the species are 
born to each alive. M 3 = I. 

(b) J of the species survive their first birthday and live into their 
second year; none of the species live longer than 2 years. In the second 
year, 3 of the species are born for each alive. M 2 = I. 

(c) j of the species survive their first birthday, of these \ survive their 
second and of these J survive their third and none of the species live 
longer than 4 years. In the fourth year 24 of the species are born for 
each alive. M 4 = I. 

8.3 100 of the species alive in the first and second year of life as a 
starting position. 

Total 
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The matrix 
tion. 



describes a species with increasing total popula- 


For the matrix ^ qJ (i) ab = \ gives a population which is steady, 

having the same values in alternate years (ii) ab> 1 the population 
increases (iii) ab< 1 the population decreases. 


8.4 P- x =_ 1 _ ( d * 

C]d 2 — C 2 d] \ ^2 ^3/ 

Now TP=T ^ and since is an eigenvector of T, then 

T (';) - A ' 0;) - (&) ” ,d s ““ y T © - ($) • 

Hence TP -(*& 

Then P _1 TP =--- ( d * (Vi X * d /\ 

c x d 2 - c 2 d x \~~ c 2 c i / \^i c 2 \d 2 ) 


. i A* 

C\d 2 — c 2 d± \ 

to- y- 


C-^2 /\ j d j f 2 0 

0 


8.5 




l_ (c x d^[ - d lC ^\ Cl d x ( a; - A n 2 )\ 

C x d 2 — C 2 d 1 1 “ ^2) — C 2^1^1 + ^1^2^27 

If A x »•= 1 and A 2 < 1 then A 2 ->0 as 00 and a finite limit for T w exists. 
If A X >1 the population increases and if A x < 1, A 2 <1 then Ai~>0, 
A^->0 as «->o 0 and T w is zero. 


8.6 (a) Steady state; eigenvalues 1 and If the start was 100 to 
each life year, the steady state is 300 and 50 in the first and 
second year of life. 

(b) eigenvalues § and - Population becomes extinct. 

(c) eigenvalues \\ and - J. Population increases. 
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8.8 C T first year 

second C /A f\ . i i j 1 
year r(j |] ^values 1 and 

Steady state of 8,000 country dwellers, 6,000 town dwellers. 

8.9 L E first lecture 

lecture E (| t) ei g envaIues 1 > " *• 

Probability of y$ late for end of term lecture. 


8.10 As matrices have probabilities as entries, the column considers 
all outcomes and then must total 1 . 


Matrices are of the form ( 1 a 
\l-a 

(a \-b\( a 1 -iW 
\1 -a b /\l-« b )~\ 

giving column totals of 1 again. 


V) and 

a 2 + 1 -a-b + ab 
a 2 -a + b -ab 


a-ab +b -b 2 \ 
\ -a-b + ab + b 2 ) 


8.12 L E 

L /° | 
E[ 1 0 

OT\ 0 | 


OT 

•§ \ Probabilities for end of term lecture: 
3 " ) L, 33 ; E, -§jj\ On Time 
0 / 


1 2 3 4 from 
1/0 f 0 0 \ 

8.13 T = to 211 0 f 0 1 
3 0 | 0 1 
4\0 0 1 0/ 

Probabilities of being in cells: 1, A; 2, fiI 3, 3 ^; 4, 5 - 4 :. 


Chapter 9 

9.1 All are 6x 2 + 4 xy + 3 y 2 = 1. 


9.2 



9.3 |)x = l; 2 ) x = l' 
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MATRICES AND THEIR APPLICATIONS 


9.4 


9.5 


pg 2 y2 

Ellipse -2 +21 = 1, x ± axis gradient 1. 


v^ 

Hyperbola -2 -22 = 1, x 1 axis gradient 


1 

2 • 


9.6 Straight lines y\ = 4, 

Xy axis gradient 2. 

9.7 Circle - note many eigenvectors obtained. 

9.8 Ellipse centre ( - 1, 1); axes 45 °, x\j 3 +y\fi = 1. 

9.9 Circle, centre (2, 1) radius 2. 

9.10 Hyperbola; centre (1, -2) axes tan -1 (|), x\j2 -jf/3 = 1. 

9.11 Parabola, axes 45 °, y\ = J2x v 



